{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [],
   "source": [
    "# Scripts to measure the orientation of a filament relative to the Galactic North \n",
    "import Orientations_Functions as OF\n",
    "import numpy as np\n",
    "from astropy.io import fits\n",
    "from aplpy import FITSFigure\n",
    "import matplotlib.pyplot as plt"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [],
   "source": [
    "List_Names = ['Fil1', 'Fil2', 'Fil4', 'Fil5', 'Fil6', 'Fil8', 'Fil10', 'G24', 'G47', 'G49']"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Fil1 angle: 92.85084909041683\n",
      "Fil2 angle: 74.08004386543406\n",
      "Fil4 angle: 95.76537653997275\n",
      "Fil5 angle: 82.0855772990879\n",
      "Fil6 angle: 86.63131057106197\n",
      "Fil8 angle: 85.62358860735564\n",
      "Fil10 angle: 94.85574563027836\n",
      "G24 angle: 107.48337412681784\n",
      "G47 angle: 69.64627733682192\n",
      "G49 angle: 98.7472400752534\n"
     ]
    }
   ],
   "source": [
    "# Measuring the angle for each filament\n",
    "Fil_Angles = []\n",
    "for i in range(0,10):\n",
    "    FileName = 'Skeletons/'+List_Names[i]+'_skeleton.fits'\n",
    "    FitsFile = fits.open(FileName)\n",
    "    # Calculating the angle from the fit parameters and the skeleton positions\n",
    "    SkeletonAngle, LineParameters, PosArray = OF.FitOrientation(FitsFile)\n",
    "    print(List_Names[i]+' angle: '+str(SkeletonAngle))\n",
    "    Fil_Angles.append(SkeletonAngle)\n",
    "    FitsFile.close()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "74.08004386543406\n"
     ]
    },
    {
     "data": {
      "text/plain": [
       "[<matplotlib.lines.Line2D at 0x1d34cc0ae00>]"
      ]
     },
     "execution_count": 4,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAA1gAAAL1CAYAAADaakUsAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjYuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/P9b71AAAACXBIWXMAAA9hAAAPYQGoP6dpAAByDElEQVR4nOzdeXiU9b3//+c9M9lXQkhYEnZl0SCCglSrVm3VWO3i75yup2p72lNrN+12errXLuecflFPN+05bbX9tp5av13ENrVaBFq3alEgWkVFESKQsIWQZbLM3L8/BiKBoETukBCej+viSub+3HPPe2ZyhXnl/bk/dxCGYYgkSZIk6bDFhroASZIkSRopDFiSJEmSFBEDliRJkiRFxIAlSZIkSRExYEmDpLW1lSuvvJLW1tahLkU6LC0tLZx99tm0tLQMdSnSYfFnWSOFnzGGNwOWNEiam5u56aabaG5uHupSpMPS0tLCihUr/FCqo54/yxop/IwxvBmwjjLJZJKvf/3rJJPJoS4FGF71DKdahqPh9voMt78kD6d6hlMtw9Fwe32s5+ioZTgabq+P9RwdtQxHw+0zxlDXY8A6ymzbto3Pf/7zbNu2bahLAYZXPcOpluFouL0+w+0vycOpnuFUy3A03F4f6zk6ahmOhtvrYz1HRy3D0XD7jDHU9SSG5FFHoG3btvHHP/6RyZMnk5eXN2iP09jYCEB9fT1NTU2D9jiHajjVM5xqAet5JdZzdNQC1vNKrOfoqAWs55VYz9FRC1jPKxmsejo6Oli/fj3nn38+5eXlB90vCMMwjOxRj2E///nPefe73z3UZUiSJEkaRD/72c9417veddBxO1gRmTx5MgAncCoFFA9tMZIkSZIi1UYLT/BI7+f+gzFgRWTvtMACiikORg1xNZIkSZIitWfe3yudDuQiF5IkSZIUEQOWJEmSJEXEgCVJkiRJETFgSZIkSVJEDFiSJEmSFBEDliRJkiRFxIAlSZIkSRExYEmSJElSRAxYkiRJkhQRA5YkSZIkRcSAJUmSJEkRMWBJkiRJUkQMWJIkSZIUEQOWJEmSJEXEgCVJkiRJETFgSZIkSVJEDFiSJEmSFBED1mFqaWmhoaGBxsbGoS5FkiRJ0iBrbGykoaGBlpaWfscTR7ieEeeSSy5hxYoVQ12GJEmSpCOgtrYWgLPOOovly5cfMG7AOkxLliyhpaWF+vr63hdbkiRJ0shUV1dHTU0NxcXF/Y4bsA5TcXExxcXFNDU1DXUpkiRJkgZZZWUlVVVVBx33HCxJkiRJiogBS5IkSZIiYsCSJEmSpIgYsCRJkiQpIgYsSZIkSYqIAUuSJEmSImLAkiRJkqSIGLAkSZIkKSIGLEmSJEmKiAFLkiRJkiJiwJIkSZKkiBiwJEmSJCkiBixJkiRJiogBS5IkSZIiYsCSJEmSpIgYsCRJkiQpIgYsSZIkSYqIAUuSJEmSImLAkiRJkqSIGLAkSZIkKSIGLEmSJEmKiAFLkiRJkiJiwJIkSZKkiBiwJEmSJCkiBixJkiRJiogBS5IkSZIiYsCSJEmSpIgYsCRJkiQpIgYsSZIkSYqIAUuSJEmSImLAkiRJkqSIGLAkSZIkKSIGrMPU0tJCQ0MDjY2NQ12KJEmSpEHW2NhIQ0MDLS0t/Y4njnA9I84ll1zCihUrhroMSZIkSUdAbW0tAGeddRbLly8/YNyAdZiWLFlCS0sL9fX1vS+2JEmSpJGprq6OmpoaiouL+x03YB2m4uJiiouLaWpqGupSJEmSJA2yyspKqqqqDjruOViSJEmSFBEDliRJkiRFxIAlSZIkSRExYEmSJElSRAxYkiRJkhQRA5YkSZIkRcSAJUmSJEkRMWBJkiRJUkQMWJIkSZIUEQOWJEmSJEXEgCVJkiRJETFgSZIkSVJEDFiSJEmSFBEDliRJkiRFxIAlSZIkSRExYEmSJElSRAxYkiRJkhQRA5YkSZIkRcSAJUmSJEkRMWBJkiRJUkQMWJIkSZIUEQOWJEmSJEXEgCVJkiRJETFgSZIkSVJEDFiSJEmSFBEDliRJkiRFxIAlSZIkSREZ1gHre9/7HpMnTyY3N5eFCxfy8MMP946tXbuW008/naqqKr72ta/1uV8ymeSqq65i9OjRFBYWcumll9LY2Nhnnw0bNnDRRReRn59PRUUFn/rUp+jp6ekdv+WWWzj77LMH9flJkiRJGlmGbcC67bbbuOaaa/jSl77Eo48+ykknncT5559PU1MTAB/+8Id597vfzR133MEdd9zBAw880Hvfq6++mjvvvJPbb7+dFStWsGnTJt761rf2jqdSKS666CK6urp44IEH+MlPfsItt9zCF7/4xSP+PCVJkiSNHMM2YF133XW8//3v54orrmD27NncdNNN5Ofn8+Mf/xiAnTt3Mn/+fObMmcP48eNpbm4GYNeuXfzoRz/iuuuu45xzzmH+/PncfPPNPPDAAzz00EMA3H333fz973/nZz/7GXPnzuXCCy/k2muv5Xvf+x5dXV1D9ZQlSZIkHeWGZcDq6upi5cqVnHfeeb3bYrEY5513Hg8++CAAX/3qVznvvPPIz88nFotx/vnnA7By5Uq6u7v73HfmzJlMnDix974PPvggNTU1VFZW9u5z/vnn09LSwhNPPHEknqIkSZKkESgx1AX0Z9u2baRSqT4BCKCyspKnnnoKgNraWrZu3UpLSwtjxozp3WfLli1kZ2dTWlp6wH23bNnSu09/x947BnD55Zdz+eWXR/m0JEmSJI1wwzJgHaqcnJw+4epIamlpoaWlpff2/otoSJIkSRp5GhsbaWho6L1dXFxMcXFx7+1hGbDKy8uJx+MHhJbGxkbGjh37svcdO3YsXV1dNDc39+li7XvfsWPH9lmRcO/43rFDcckll7BixYpD2leSJEnSyFBbW9vn9llnncXy5ct7bw/Lc7Cys7OZP38+S5cu7d2WTqdZunQpixYtetn7zp8/n6ysrD73Xbt2LRs2bOi976JFi6ivr+9dkRDgnnvuobi4mNmzZx9SjUuWLGHjxo29/+rq6gbyFCVJkiQdherq6vrkgCVLlvQZH5YdLIBrrrmGyy67jFNOOYUFCxZwww030NbWxhVXXPGy9yspKeF973sf11xzDWVlZRQXF/ORj3yERYsWcdpppwHwhje8gdmzZ/NP//RP/Od//idbtmzh85//PFdddRU5OTmHVN/+rcB9w5okSZKkkamyspKqqqqDjg/bgPW2t72NrVu38sUvfpEtW7Ywd+5c7rrrrgMWp+jP9ddfTywW49JLL6Wzs5Pzzz+f73//+73j8Xic3/3ud1x55ZUsWrSIgoICLrvsMr761a8O5lOSJEmSNMIFYRiGQ13ESPDoo48yf/58FnAuxcGooS5HkiRJUoRawp08zFJWrlzJvHnzDrrfsDwHS5IkSZKORgYsSZIkSYqIAUuSJEmSImLAkiRJkqSIGLAkSZIkKSIGLEmSJEmKiAFLkiRJkiJiwJIkSZKkiBiwJEmSJCkiBixJkiRJiogBS5IkSZIiYsCSJEmSpIgYsCRJkiQpIgYsSZIkSYqIAUuSJEmSImLAkiRJkqSIGLAkSZIkKSIGLEmSJEmKiAFLkiRJkiJiwJIkSZKkiBiwJEmSJCkiBixJkiRJiogBS5IkSZIiYsCSJEmSpIgYsCRJkiQpIgYsSZIkSYqIAUuSJEmSImLAkiRJkqSIGLAkSZIkKSIGrMPU0tJCQ0MDjY2NQ12KJEmSpEHW2NhIQ0MDLS0t/Y4njnA9I84ll1zCihUrhroMSZIkSUdAbW0tAGeddRbLly8/YNyAdZiWLFlCS0sL9fX1vS+2JEmSpJGprq6OmpoaiouL+x03YB2m4uJiiouLaWpqGupSJEmSJA2yyspKqqqqDjruOViSJEmSFBEDliRJkiRFxIAlSZIkSRExYEmSJElSRAxYkiRJkhQRA5YkSZIkRcSAJUmSJEkRMWBJkiRJUkQMWJIkSZIUEQOWJEmSJEXEgCVJkiRJETFgSZIkSVJEDFiSJEmSFBEDliRJkiRFxIAlSZIkSRExYEmSJElSRAxYkiRJkhQRA5YkSZIkRcSAJUmSJEkRMWBJkiRJUkQMWJIkSZIUEQOWJEmSJEXEgCVJkiRJETFgSZIkSVJEDFiSJEmSFBEDliRJkiRFxIAlSZIkSRExYEmSJElSRAxYkiRJkhQRA5YkSZIkRcSAJUmSJEkRMWBJkiRJUkQMWIeppaWFhoYGGhsbh7oUSZIkSYOssbGRhoYGWlpa+h1PHOF6RpxLLrmEFStWDHUZkiRJko6A2tpaAM466yyWL19+wLgB6zAtWbKElpYW6uvre19sSZIkSSNTXV0dNTU1FBcX9ztuwDpMxcXFFBcX09TUNNSlSJIkSRpklZWVVFVVHXTcc7AkSZIkKSIGLEmSJEmKiAFLkiRJkiJiwJIkSZKkiBiwJEmSJCkiBixJkiRJiogBS5IkSZIiYsCSJEmSpIgYsCRJkiQpIgYsSZIkSYqIAUuSJEmSImLAkiRJkqSIGLAkSZIkKSIGLEmSJEmKiAFLkiRJkiJiwJIkSZKkiBiwJEmSJCkiBixJkiRJiogBS5IkSZIiYsCSJEmSpIgYsCRJkiQpIgYsSZIkSYqIAUuSJEmSImLAkiRJkqSIGLAkSZIkKSIGLEmSJEmKiAFLkiRJkiJiwJIkSZKkiBiwJEmSJCkiBixJkiRJiogBS5IkSZIiYsCSJEmSpIgYsA5TS0sLDQ0NNDY2DnUpkiRJkgZZY2MjDQ0NtLS09DueOML1jDiXXHIJK1asGOoyJEmSJB0BtbW1AJx11lksX778gHED1mFasmQJLS0t1NfX977YkiRJkkamuro6ampqKC4u7nfcgHWYiouLKS4upqmpaahLkSRJkjTIKisrqaqqOui452BJkiRJUkQMWJIkSZIUEQOWJEmSJEXEgCVJkiRJETFgSZIkSVJEDFiSJEmSFBEDliRJkiRFxIAlSZIkSRExYEmSJElSRAxYkiRJkhQRA5YkSZIkRcSAJUmSJEkRMWBJkiRJUkQMWJIkSZIUEQOWJEmSJEXEgCVJkiRJETFgSZIkSVJEDFiSJEmSFBEDliRJkiRFxIAlSZIkSRExYEmSJElSRAxYkiRJkhQRA5YkSZIkRcSAJUmSJEkRMWBJkiRJUkQMWJIkSZIUEQOWJEmSJEXEgCVJkiRJETFgSZIkSVJEDFiSJEmSFBEDliRJkiRFxIAlSZIkSRExYB2mZDJJQ0MDjY2NQ12KJEmSpEHW2NhIQ0MDLS0t/Y4bsA7T4sWLqa6upra2dqhLkSRJkjTIamtrqa6u5pJLLul3PAjDMDzCNY0oyWSSbdu2UV9fT21tLQs4l+Jg1FCXJUmSJClCLeFOHmYpdXV11NTUUFxcTHFx8QH7JYagthElNzeXqqoqmpqahroUSZIkSYOssrKSqqqqg447RVCSJEmSImLAkiRJkqSIGLAkSZIkKSIGLEmSJEmKiAFLkiRJkiJiwJIkSZKkiBiwJEmSJCkiBixJkiRJiogBS5IkSZIiYsCSJEmSpIgYsCRJkiQpIgYsSZIkSYqIAUuSJEmSImLAkiRJkqSIGLAkSZIkKSIGLEmSJEmKiAFLkiRJkiJiwJIkSZKkiBiwJEmSJCkiBixJkiRJiogBS5IkSZIiYsCSJEmSpIgYsCRJkiQpIgYsSZIkSYqIAUuSJEmSImLAkiRJkqSIGLAkSZIkKSLDMmCFYcgXv/hFxo0bR15eHueddx7PPPNMn30efPBB5s6dy+TJk/nRj37UZ+zrX/86r3nNa8jPz6e0tLTfxwiC4IB/v/jFL3rHb7nlFs4+++yon5okSZKkEWxYBqz//M//5Nvf/jY33XQTf/3rXykoKOD8888nmUz27vO+972PL3zhC9x6661885vfZOPGjb1jXV1d/MM//ANXXnnlyz7OzTffzObNm3v/vfnNbx6spyRJkiTpGJAY6gL2F4YhN9xwA5///Od505veBMBPf/pTKisr+e1vf8vb3/52ANra2pg3bx4VFRWMGjWK3bt39x7jK1/5CpDpQr2c0tJSxo4dOzhPRJIkSdIxZ9h1sJ5//nm2bNnCeeed17utpKSEhQsX8uCDD/Zu++IXv8isWbMoKSnhtNNOY/bs2QN+rKuuuory8nIWLFjAj3/8Y8IwjOQ5SJIkSTo2DbsO1pYtWwCorKzss72ysrJ3DDJTBN/+9rfT1dXFqFGjBvw4X/3qVznnnHPIz8/n7rvv5kMf+hCtra189KMfBeDyyy/n8ssvf/VPRJIkSdIxZ8gD1s9//nP+5V/+pff273//+0O+b0FBAQUFBa/qcb/whS/0fn/yySfT1tbGt771rd6A9UpaW1tpbm7uvd3Y2Piq6pAkSZJ09GhsbKShoaH3dmlpKYWFhb23hzxgXXLJJSxcuLD3dmdnJ5ApfNy4cb3bGxsbmTt37qDVsXDhQq699lo6OzvJycl5xf0/9alPcdNNNw1aPZIkSZKGn9ra2j63P/jBD3LjjTf23h7ygFVUVERRUVHv7TAMGTt2LEuXLu0NVC0tLfz1r399xVUBD8eqVasYNWrUIYUrgG9961t87nOf671dX19/wIstSZIkaWSpq6ujpqam9/b+l4Ua8oC1vyAI+PjHP87XvvY1jjvuOKZMmcIXvvAFxo8ff8jLqG/YsIEdO3awYcMGUqkUq1atAmD69OkUFhZy55130tjYyGmnnUZubi733HMP3/jGN/jkJz95yHUWFhb2aQU2NTUN5GlKkiRJOgpVVlZSVVV10PFhF7AAPv3pT9PW1sYHPvABmpubOeOMM7jrrrvIzc09pPt/8Ytf5Cc/+Unv7ZNPPhmAZcuWcfbZZ5OVlcX3vvc9rr76asIwZPr06Vx33XW8//3vH5TnI0mSJOnYEISuTR6JRx99lPnz57OAcykOBr6qoSRJkqThqyXcycMsZeXKlcybN++g+w2762BJkiRJ0tHKgCVJkiRJETFgSZIkSVJEDFiSJEmSFBEDliRJkiRFxIAlSZIkSRExYEmSJElSRAxYkiRJkhQRA5YkSZIkRcSAJUmSJEkRMWBJkiRJUkQMWJIkSZIUkVcVsJqbm/nhD3/IZz/7WXbs2AHAo48+yosvvhhpcZIkSZJ0NEkM9A5r1qzhvPPOo6SkhPXr1/P+97+fsrIyfv3rX7NhwwZ++tOfDkadkiRJkjTsDbiDdc0113D55ZfzzDPPkJub27u9traWP//5z5EWJ0mSJElHkwEHrEceeYR/+Zd/OWD7hAkT2LJlSyRFSZIkSdLRaMABKycnh5aWlgO2P/3004wZMyaSoiRJkiTpaDTggHXJJZfw1a9+le7ubgCCIGDDhg185jOf4dJLL428QEmSJEk6Wgw4YC1evJjW1lYqKiro6OjgrLPOYvr06RQVFfH1r399MGqUJEmSpKPCgFcRLCkp4Z577uG+++5jzZo1tLa2Mm/ePM4777zBqE+SJEmSjhoDDlh7nXHGGZxxxhlR1iJJkiRJR7VDCljf/va3D/mAH/3oR191MZIkSZJ0NDukgHX99df3ub1161ba29spLS0FoLm5mfz8fCoqKgxYkiRJko5Zh7TIxfPPP9/77+tf/zpz587lySefZMeOHezYsYMnn3ySefPmce211w52vZIkSZI0bA14FcEvfOELfOc732HGjBm922bMmMH111/P5z//+UiLkyRJkqSjyYAD1ubNm+np6TlgeyqVorGxMZKiJEmSJOloNOCAde655/Iv//IvPProo73bVq5cyZVXXulS7ZIkSZKOaQMOWD/+8Y8ZO3Ysp5xyCjk5OeTk5LBgwQIqKyv54Q9/OBg1SpIkSdJRYcDXwRozZgx1dXU8/fTTPPXUUwDMnDmT448/PvLiJEmSJOlo8qovNHz88ccbqoBkMsm2bds8/0wj1vagkdFh5VCXIUmSNCw0NjbS0NBAeXk5ubm5B4wPOGC9973vfdnxH//4xwM95FFt8eLFrp6oEas77GJt0SoWtJxDIsga6nIkSZKGXG1tLQBf+9rX+NznPnfA+IAD1s6dO/vc7u7u5vHHH6e5uZlzzjnnVZZ59PrEJz7BZZddRn19fe+LLY0UW7M3E56aZutfNpHdk2snS5IkHfPq6uqoqamhvLy83/EBB6zf/OY3B2xLp9NceeWVTJs2beAVHuVyc3OpqqqiqalpqEuRIre9dAvxqVk0rd5EW1eLnSxJknTMq6yspKqq6qDjr/ocrH3FYjGuueYazj77bD796U9HcUhJR1hnmGR90VO8mPM8AFnTszn5hDTdJ8ETy7pJzmhj+UNLIAXVqWlM3j2TnODAeceSJEnHskgCFsC6dev6vQCxpKNDNjnk9xRROK2Y9CkhQSyge047AD0VXQRzYsSaYoQvpMlpySebnCGuWJIkafgZcMC65ppr+twOw5DNmzfz+9//nssuuyyywiQdWUEQUJ2cRsGqIp7Y8giddLByCQQnBORXQDoXkrvThKemefahep4rfIIJnVPsZEmSJO1jwAHrscce63M7FosxZswYFi9e/IorDEoa/kYxhom7prNu7t8Jz8p0ssL8NLE2iJfFSc2B+I44WVuzyX+xyE6WJEnSPgYcsJYtWzYYdUgaJoIgIK+7iKxYNp3/L0nYmqY9AGogKEwRhAHpZJqO3e08nbeaZwrW2MmSJEnaIzbQO5xzzjk0NzcfsL2lpeWYXKZdGom2lrxIYkyC/OZCYlNicDmwEGKnpckL0kw4OXM7+5QcCquKye+xkyVJkgSvImAtX76crq6uA7Ynk0n+8pe/RFKUpKGTDlPsTjdT+aeJTG8+key8PGJL47AF0mlIJWDH+jRUhaSfSjN79SlUJ6cRBMFQly5JkjTkDnmK4Jo1a3q///vf/86WLVt6b6dSKe666y4mTJgQbXWSjrg0aWY1z6eEMp4oe4RYTkD24zkk17cTxiBZE0IKgkKoOjnFY+Pvo/OpJLFEzKmCkiTpmHfIAWvu3LkEQUAQBP1OBczLy+M73/lOpMVJOvISQRYllL3UyfrzRDrpYMuCFwgXAN0hYRhCHmydkSKrtICuTZ0UlBSR/6xTBSVJ0rHtkAPW888/TxiGTJ06lYcffpgxY8b0jmVnZ1NRUUE8Hh+UIiUdeXs7WblhPivHLCddn4b1wHlAQWafZDqkZ8sumA27/9rM9lgjFUwgBztYkiTp2HTIAWvSpEkApNPpQStG0vCxt5OVooeedA85k3JJzM6i49l2WATxIEVPHrApIKiF2LNxdm3dTmzgp3ZKkiSNGIcUsJYsWcKFF15IVlYWS5Ysedl9L7nkkkgKkzQ8hISctOM1xHfEWfPCQ6SnpuDmzIIXzIFgbEhuT0CKBF1nJFnx0J3EAs/HkiRJx6ZDClhvfvOb2bJlCxUVFbz5zW8+6H5BEJBKpaKqTdIwsLeTBVAYltCe2wrTQjgtgFTmXKw8AloqUgRzYl6EWJIkHdMOKWDtOy3QKYLSsSkdpmhhB1ntWcQq4nTenoROCAPYcUKaWCFehFiSJB3zBnyyxE9/+lM6OzsP2N7V1cVPf/rTSIqSNPykSVOQLiZnYi5FfyulqmkqTAXeCyyEnNeG5AVpsiqyvAixJEk6Zg04YF1xxRXs2rXrgO27d+/miiuuiKQoScNPjBidQQcV91Rz0tbXUNZdSW5uAdwTQBMQZi5C3N3c7UWIJUnSMeuQVxHcKwzDfj8sNTQ0UFJSEklRkoaffS9ATABbS14klhOQ+3geyfXtdMRCqIFYQcqLEEuSpGPWIQesk08+ufdCw+eeey6JxEt3TaVSPP/881xwwQWDUqSkobfvYhcvdxHiNCGBFyGWJEnHqEMOWHtXD1y1ahXnn38+hYWFvWPZ2dlMnjyZSy+9NPICJQ0/L3cR4lgxhBx4EeL2xG666PQixJIkaUQ75ID1pS99CYDJkyfztre9jdxcPyRJx6qXuwhx7rlpunvC3osQ8/qQxPNZ5G+wgyVJkka+AZ+Dddlllw1GHZKOQv1dhLj9e2RaWHOAihA6ICc3l/UnPcUzz64B8HwsSZI0Yg04YKVSKa6//np++ctfsmHDBrq6uvqM79ixI7LiJA1vr3QR4iAOQQxSp6UIXoyRf2IRsb8F5K+1myVJkkamAS/T/pWvfIXrrruOt73tbezatYtrrrmGt771rcRiMb785S8PQomShrt9L0KcO6aA2B9iBN0x6IYwG7p7uuhc1UHi91ku3S5Jkka0AQesn//85/zP//wPn/jEJ0gkErzjHe/ghz/8IV/84hd56KGHBqNGScNcmjR5PQWEqZDkfe2Em9OEj6QJsiEWg1RTirAmTXPnVh4efS9rC1fRGSaHumxJkqTIDThgbdmyhZqaGgAKCwt7Lzr8xje+kd///vfRVifpqBAjRlesk5LnyymeUErOafkEs2PEsiBdD+wMYQrEquIUVhWT3+MUQUmSNDINOGBVVVWxefNmAKZNm8bdd98NwCOPPEJOjh+YpGNRmjSzd53CyS2nM31VDV2rkoR/SNPzQwibIUgDBZDOSZHMbeeZgjUsK/+tnSxJkjTiDHiRi7e85S0sXbqUhQsX8pGPfIR3v/vd/OhHP2LDhg1cffXVg1GjpGFu38UuSigj3pkgnJ6G1wUEsYCgK0UqG1gAOd15ZJ+Ii11IkqQRacAB69///d97v3/b297GpEmTeOCBBzjuuOO4+OKLIy1O0tEnSQf5PYXsWrcdmkLC8yAoAeJAVsDurmb4E7AroMSLD0uSpBFmwFME93faaadxzTXXsHDhQr7xjW9EUZOko1g+hYSJNLGaGLGz4sSejUEuEAAFIfwtoOCMYoonlXouliRJGnEOO2DttXnzZr7whS9EdThJR6kk7bQHbaRXp0nfnSI9OQ0p4AmgFdgW0vanFlpeaKZ9TwdLkiRppIgsYEkSQIIsJrRNpnjiKOJVCUhB8CDkjAHuBeZDUBojb3S+HSxJkjTiGLAkRSoryOa45BxmrppH+vk0PAfZZ0BsRwBlwBQId6cJm0LGdUz0gsOSJGlEGfAiF+ormUyybds2Ghsbh7oUaViJkyArnUPXsx0knwZOgFgIYR7klENybjvLH1pCLIhR1lHJrLZ55AQudiFJkoa3xsZGGhoaKC8vJzf3wM8uhxywrrnmmpcd37p168CrGwEWL17M5z//+aEuQxp2cskjTQ9MCwjOCAi6A3JK03S3hHQVQzAnRnxHnLzd+ezauJ2YDXVJknQUqK2tBeBrX/san/vc5w4YD8IwDA/lQK973esO6QGXLVs2gPKOfns7WPX19dTW1rKAcykORg11WdKQ6wm7qS/9K201u+h6rpN0Kp0ZqAFikLsoIOc+aN0UJ1XVAysDqrumMXn3TDtZkiRp2GkJd/IwS6mrq6OmpubwO1jHWnA6VLm5uVRVVdHU1DTUpUjDSowYnUEH4x6dQjwd57m5fyc8K4R4SEhIPAmTL87iidtCgvkx4i/ECV8IXfRCkiQNa5WVlVRVVR103HOwJA2KNGlmNc+nhDKa4pvIzsul864kYV4a1kMbIY/P6SLMixFriZGdlcuLJc/xYv7zBLGACZ1T7GZJkqSjjgFL0qBIBFmUUAbA1pIXSZTEyV9dzu7iXXTPScIpkCYgPx0yfWyMp55NEU6FnGQOOU255K91CXdJknT08axySYMqHabYnW6m4p5qqnZNJbsmm7zCArgjIGwJSSWgtT0k1dwNVSHpp9LMXn0K1clpLuEuSZKOOnawJA2qfacKPlH2CEEc8pcXUtBVwrbbNpGMh6w7oZvsYggKIF4W57Hx99H9bBeAUwUlSdJRxYAlaVDtnSq4bydrcusMnip4jNgJMdIL09AFqQDCbOh+TSf5PcVknZhD7G+BUwUlSdJRZcAB6+abb6awsJB/+Id/6LP99ttvp729ncsuuyyy4iSNHPt2sjpJ0py/lXR9ZsELzoOwCIhDKjvF7hd2whPAroCSxG666CQHO1iSJGn4G/A5WN/85jcpLy8/YHtFRQXf+MY3IilK0siz76IXCRL0pHvImZhLwelFxJ6NExQCARAHVkL+giKKJ5WS32MHS5IkHT0G3MHasGEDU6ZMOWD7pEmT2LBhQyRFSRrZknRQ2lVO0/Mv0rklCedD0A3BsxCmgEXQvnw3dNvBkiRJR5cBd7AqKipYs2bNAdtXr17N6NGjIylK0siWTyGjesZQPKmUeFUCUhA8CNnlwNPADAhGBxRPKiXs9uLDkiTp6DHgDtY73vEOPvrRj1JUVMSZZ54JwIoVK/jYxz7G29/+9sgLlDTyBEFAdXIaeasKWFVxH7wA6Rh0BkA10AVhfkh70EZL3k5ezH+equ6priYoSZKGvQEHrGuvvZb169dz7rnnkkhk7p5Op3nPe97jOViSBqSEMvJ6Cuk4uRVOhCAnRiydJp0NRecGdK6Ok5oYI2dNLmFDppO1PWhkdFg51KVLkiT1a8ABKzs7m9tuu41rr72W1atXk5eXR01NDZMmTRqM+iSNYHHihKSJNyYIW0LCuSHZ1QHdLSHJ7pDUzm6CaQF5fy5kR3YTyWQHawtXsaDlHBJB1lCXL0mSdIBXfR2s448/nuOPPz7KWiQdY9KkKUgX03lcO1krctjdtIuOzq7MaoInQEFxiukzslhb2kZyRhv3P1gHp8Cahx5kfOsUxobVAHa1JEnSsHFIAeuaa67h2muvpaCggGuuueZl973uuusiKUzSyBcjRmfQQcU91RR0FPHUcY8RK46RPiEN+ZCKQ2t7SE9FF8GcGDwfwsyArtVJnuPvlLeMJSRkbZFdLUmSNDwcUsB67LHH6O7u7v1ekqKw78WHnyh7hJzxOcSXF9Je30p3qpNkPGRdmPndQyHEyiGWE9I6qgVOheUP3gEEJObF2bTiBSZ2Tx/KpyNJknRoAWvZsmX9fi9Jh2PvxYfTYYrd6WYq7qlmcusMNuY9w7MnPUE4Ng3dEMyDIICcrIAggLAqJHUSxLbHiO2IE5sSZ9fq7WzfVsTosLLPlMFDmT7oFENJkhSVAV8H673vfS+7d+8+YHtbWxvvfe97IylK0rFlbydratssYkGMvO4icvLziG1KQAzC2yH9P9DxnZD2/w1J1QM/hnRWmtSYHroSSRrHbeSxMX/hTyW/4rGSv/BkwUrawt2sLVxFz94uWD+6w65X3EeSJOlQDThg/eQnP6Gjo+OA7R0dHfz0pz+NpChJx5a9nay9tpa8SKIkzqjnysl/tphgfADHA+8HZgDzIXgf5L0+IO8CyB8XwMQATofEtATBohjpMM2u7O2Ep6bZmr2J7UEjwAFft2ZvPmAfSZKkV+uQA1ZLSwu7du0iDEN2795NS0tL77+dO3dSV1dHRUXFYNYq6Riw73TBCc1TCI6DnOw8KASWA2OAKuDPkGoOCVOQaIP49jjBcTHSO1LkTcgnXZBme2kj8alZNBVvYm3hKjrC9j5fe8Jutpdu6bOPnSxJknQ4DnmZ9tLSUoIgIAiCfpdnD4KAr3zlK5EWJ+nYs//CF0EcwvqQvAUFdD3fSTqdJpyRJkxD10+AREhHCKR6YBeE7dD+u1baZ7ZCGsiG1tHNMB3uv78OJsD9TXUQg/tOqOOU2TG6T4InlnWTnNHG8oeWQAqqU9OYvHsmOUHu0L4gkiTpqHLIAWvZsmWEYcg555zDr371K8rKXprOk52dzaRJkxg/fvygFCnp2LH/whdj7p7A1uxNVNxTRX6yiLULH6M73glbApgbZjpbSeAU4CkIksDJEMSBNMRj0F0RwJwQngMmAgUQnBODFuienXncvUvBx5pihC+kyWnJJ5ucIXoVJEnS0eqQA9ZZZ50FwPPPP8/EiRMJgmDQipKkvZ2sgrCI8vZxLy3lTg4FdxYRS8XZcVwTPBXC8QH8KYSFEO4EfrnnawzSPUA8hKf3HLgKSEL46zQ9LWlW/gY4ASiCeBxSu4FT0zz7UD3PFT5BWUcls9rm0Rrb5UqDkiTpFR1ywNrr3nvvpbCwkH/4h3/os/3222+nvb2dyy67LLLiJB27ehe+CDhgKfdJrcfxyKhlTPnzTFJhiobsdeTNLaC9fjfhaSFsAiqA0UAnkLfn6zyIxSEMIXwM2AWcCQSZE1LjHRCUxUnNgfiOOHm789m1cTtpUqwt9GLGkiTplQ14FcFvfvOblJeXH7C9oqKCb3zjG5EUJUn723cp95CQWc3zmdZ+AqVd5WRPzqWnqQcqA+gKYDdQDjSRuUDxdmAMBH+GdAsEISRKgGLgNuC/M8vAd62GVGGKIAwIe0La0210L+zk/rI/0LWwg03ZLwzhKyBJko4GAw5YGzZsYMqUKQdsnzRpEhs2bIikKEna375Lue/7/daSF0nE41AfkpvMI/u2bGJtcciCWEuc+FNZBNMS8CKE84HVQE+c1PoAKgKCtoDYlBhcDiyE2Glp8oI0J747RlAaEMyPEZQFxKYk2FW6/YBl3iVJkvY14CmCFRUVrFmzhsmTJ/fZvnr1akaPHh1VXZL0ivpbCKO6dRoPtv+RcEucUR1j2Hr+ZrI35ZIq74GNwDhIt6fgBYgHCeITE4T5IeEfugnz0qTWQVssZHXYlXmQfIiNidFd0EnjuI00Bhthz9CE7ilMbT3BlQYlSVKvAQesd7zjHXz0ox+lqKiIM888E4AVK1bwsY99jLe//e2RFyhJB9PfQhg9dHPSjtMB2Fj2LLlhHun6NHkLCkiubCc8G/gVZI3NJvZcjFRlN1TBqFVj2F28i+65SSgAuiA4GXJzAoIgBKD9uACmhwTrY1AJ6RVpsslhe9DoAhiSJAl4FVMEr732WhYuXMi5555LXl4eeXl5vOENb+Ccc87xHCxJR9TeqYL9TR8sCktoTe9i7N0TyWrPZuzdE8nfXcSoO8vJbymk6smpxHcnCNbHiS2LM37XFLJrsskrKoBnAhhF78WMAcIUBJuB6QHBDnovZtxDtxcoliRJvQbcwcrOzua2227j2muvZfXq1eTl5VFTU8OkSZMGoz5JelX2727t2+XaxQ4KwmJGMQaCgDhxNpY9SxCHwr+UMLNpHqvTD5AuS9H1E+hKhJDJWQT5kDcNuka10jSplcYxG2ECLF+5hPLkOGa1zXPKoCRJx7ABB6y9jj/+eI4//vgoa5GkyOy/zHvvV166XUYFkDmXqzW9i4p7qpnSOpOtic3kTsunfXNr5gLFJbx0MeMY9MQgKx6QdQm03w3MD4k/F2fX1u3EBj4xQJIkjSCvKmA1NDSwZMkSNmzYQFdXV5+x6667LpLCJOlI2dvt2hu8tpa8SCwRkLM5l3Rpiu413TAvhDoI86Hz2T3trBMgtxiyw4AkCbrOSLLioTuJBTEmdE5h8u6ZdrMkSTrGDDhgLV26lEsuuYSpU6fy1FNPceKJJ7J+/XrCMGTevHmDUaMkDap9z+E6YGXCh6oY3zqZR/LupfOUDngMmAecBDHixMM0U8Zm8URFimBOjPiOOFlbs8l/sYhscob0eUmSpCNvwAHrs5/9LJ/85Cf5yle+QlFREb/61a+oqKjgXe96FxdccMFg1DisJZNJtm3bRmOj18SRRoL+ViZsSmwiMTWL9JaQntldhGEa/hfSHSnaAnh8ThdhYZwgDEgn03TsbufpvNU8U7DGTpYkSSNMY2MjDQ0NlJeXk5t74P/vAz5Z4Mknn+Q973kPAIlEgo6ODgoLC/nqV7/Kf/zHfxx+xUeZxYsXU11dTW1t7VCXIikC/a1MuLXkRYI4xOoDcpN5ZN2XTbwgC6YC74X0KQG5Z2YuUJxVkUVwekD2KTkUVhUTdod2siRJGkFqa2uprq5m8eLF/Y4POGAVFBT0nnc1btw41q1b1zu2bdu2V1nm0esTn/gEGzdupK6ubqhLkTQI9p0ymNWezbh7JnNGw0WUrx+buV7WPQFhU0iYhlQCenZ3w3EQrk1z/KqT2JHdRIqeoX4akiQpInV1dWzcuJFPfOIT/Y4PeIrgaaedxn333cesWbOora3lE5/4BPX19fz617/mtNNOO+yCjza5ublUVVXR1NQ01KVIGgT9TRkkgFhOjNxEHuF66FzfQUcshBpIjEqRHwuYeTaseuR+eqq6Wb5yCdVd05wqKEnSCFBZWUlVVdVBxwccsK677jpaW1sB+MpXvkJrayu33XYbxx13nCsIShpxDljunZe6WuP+PIVOOtiy4AXCBUB3SCo/JBVA88w0PAUcB4n1CfI3uuiFJEnHggEHrKlTp/Z+X1BQwE033RRpQZI03O27rPsTZY+QE8+j+9YuwsqQVEMPySBk3QndZBVBrADi+QnWTXmcZ5rWkO5JU52ymyVJ0kjlFTElaYD2drV6O1kPTWbRxvPJ2p5NbE4MLgcWQroGcooCeuZ3k2pNkXteAbHsgJzufLtZkiSNUIfUwRo1ahRBEBzSAXfs2HFYBUnS0WLfTtaL8fX0FHWTfiwNzwCvh3QxpNOQakrBKSHty1pgLux4aAvjkhPJwQ6WJEkjzSEFrBtuuGGQy5Cko8++S7nvLGkiVhAjMTmLxIxsks+2w2tDelIhtAZwMgRPBjAzoLu+i92tzeSEY4f4GUiSpKgdUsC67LLLBrsOSTpqpcMUreEuquqnU9JaRv0Lf4UwJPw7pEKAEBIQloUEcdhd0cyq9H3QA7Eg7sWIJUkaQQa8yMW+kslk7zWx9iouLj6sgiTpaLPvVMGQkLxUPrtP6oJcIATmAXGIAznZ0DMbuqZDvDVBQWMR4VNejFiSpJFiwItctLW18eEPf5iKigoKCgoYNWpUn3+SdKzZd6pgSJrusIt4U4LY1jgUQuw+YFfmfKx0D6Q3A9MD4k/FvRixJEkjzIAD1qc//WnuvfdebrzxRnJycvjhD3/IV77yFcaPH89Pf/rTwahRko4aadIUpIvJOy6fUc+PIW91EWEygF9A+CNI3heSzo5REIOas2HV7PtpP203y0cvYVX+A3SGSbYHjb3H2/v9vtskSdLwNeCAdeedd/L973+fSy+9lEQiwWtf+1o+//nP841vfIOf//zng1GjJB01YsToDDqouKeaCc1TCGZAvCcBJwLvBU6GcD6kEnsuRlwKwfwY8VFxdiW2kybF2sJV9ITddIddrC1cRUfY3rtNkiQNbwMOWDt27Oi92HBxcXHvsuxnnHEGf/7zn6OtTpKOMnvPx5raNoutJZuIJQLiWxJkB7nwM+A+CG9Ok/xhyLql3WSX9VAUhsRJ0H1GJ/eX/YH203azpvhBNmdvIDw1zfP5TxKemmZr9qahfnqSJOkVDDhgTZ06leeffx6AmTNn8stf/hLIdLZKS0sjLU6Sjjb7X4R4zN0TyGrPpuqhaZz+Qi3xpiw4id6LEQcnwpRxWaQrUgRzYgRjA4KZMboLuthVup341Cy2lzcSn5rFttItThmUJGmYG/AqgldccQWrV6/mrLPO4l//9V+5+OKL+e53v0t3dzfXXXfdYNQoSUedvZ2sgrCI8vZxlFBGU2ITuSfk0ZWK031rJ7SHtMVDVoddkAJ2AU3AL0N2h83sDprheaAaOgvbaR3XTGOwEfYs3jqhewpTW09weXdJkoaRAQesq6++uvf78847j6eeeoqVK1cyffp05syZE2lxknS06l1ZMKB3hcGtJS8SxKHkwVHktubTMHcdjCUTmE4GngQmAzkvbQviEKQhEYOu4wKYHhKsj0ElpFekySaH7UEjo8PKIXqmkiRpXwOeIri/SZMm8da3vtVwJUkvY++UwYp7qjlp62so664kN68AXgwgH2J/ASZC1hSIbQdGQ7ACaIbsREA8CAj2LO8e7IC8CfmkC9L00O0CGJIkDSOH3MHq6Ohg6dKlvPGNbwTgs5/9LJ2dnb3j8Xica6+9ltxcp6pI0v72vRgxQaablSiJM3pNJW2tu+k+rRMe76FnQUAwIyB4HMIwDT+DZCLMXLAYCPIhbxp0jWqlaVIrjWM2wgRYvnIJ5clxzGqbR2tsV29Ha293yy6XJElHxiF3sH7yk5/wgx/8oPf2d7/7XR544AEee+wxHnvsMX72s59x4403DkqRknS02/dixPt2syY0TyE+PU7Wpmwoh7ArJL06TXhcOnPHOcAiYB7wAeA10PNayCoOyLskgLIA5uMy75IkDROHHLB+/vOf84EPfKDPtltvvZVly5axbNkyvvWtb/WuKChJOrj9l3IP4hDWh+QlC0j8IotYWxw6INgeEGTFCNYEUAjUQXgvdN4U0v69kPYVIbnFIcVhcMjLvPe3CuH+2wY6JkmSXnLIAevZZ5+lpqam93Zubi6x2Et3X7BgAX//+9+jrU6SRqCDLeU+7p7JnL7xQhKbEySWZpHYnk3WyiwWrbuAnCfzoAZoIdPNugJiC+LE5wSHvMx7U/GmAzpa+3e5BjpmV0ySpL4O+Rys5ubmPudcbd26tc94Op3uMy5Jenn9LeXeQzcn7TidVNgDQUCcOC2JHSSmZpHeEtIzuytzbtb/Qro1RVucQ1/mfXQzHAf3P1AHp8Cahx6kPDk+0+V68Mk+X7f+ZRMhvOLY1r9sYlz3pKF7ESVJGmYOOWBVVVXx+OOPM2PGjH7H16xZQ1VVVWSFSdJI199S7vtu26uhZB1BHGL1AbkL8uh5pJt0aUhqbPfAlnkfH8BJITxHpstV38Wu3D1drmcayZqa3ft1W/0WgFcc21a/hextuS6kIUnSHoccsGpra/niF7/IRRdddMBKgR0dHXzlK1/hoosuirxASTqW7TuNcGv2JiruqWJS6/H8vflvNC7amFnmfVJI7C+QngdZ3ZCqh3Ql8EsId0IYg64QSIWZDlcThL9Ms7unmd3xZlgCpKDz/nZoh84722lNN0OcvmOd0Pm7dlpnNkMayOKlix93BhCEXvxYknTMO+SA9W//9m/88pe/ZMaMGXz4wx/m+OOPB2Dt2rV897vfpaenh3/7t38btEIl6VjU3zRCAojlxMgrLCDvsULaW1sPXOb9aQir01ABjAE6ObDL1clL3a65wFPAlJcZS9KnKxaPQffEAKaFJF5IkKpMkV6RZnesmZxw7BF8lSRJGj4OOWBVVlbywAMPcOWVV/Kv//qvhGHmoixBEPD617+e73//+1RWOjVEkqLU3zTCvV2tcfdMpjBZzLrXPAGbskmV9xB2hYSrQ5gVZBbE6AQagSoI/gLhvl2uyRA8D+GUPWPzM/sH9Xu29Te2ApgHWSWQkwxIb4+Tfm2a9GMp8ubl05PXw9Px1ZS2jCYRZA3VyyZJ0pA55IAFMGXKFO666y527NjBs88+C8D06dMpKysblOIkSQfa96LFT5Q98tIy7wsK6P5FF+mSNOn2FMHjAbE5cWKPxehJ9xDOTsPj9Ha5Ys8EMDMkeCYgnBfC6jSxhfHebfuPBSfECLszFz/uTEBnGAI9BAWQNy2ga1QrW8tbYTosf/AOYkGcso5KZrXNc8qgJOmYMaCAtVdZWRkLFiyIuhZJ0iHob5n3vednVbdO48H2P5LeHCMgRrAcjm+ey9qFq+je3Nl7MePeLte2EMYG0BDCOEh3peEgY2FWGprJTBksJDNl8FQgBj0xyIoHpKtCUidBbHuMgt2F7Nq4ndihXxFEkqSj3qsKWJKkoXeoy7xvLHuWHHJgNSQWxOn+RfcBXa7w4TTps0OC2wOColj/Y/8vICgKSGelMudkTQR+AWErdBLCiUAcchKQWxDSuquN1MIeVqy8k+quaUzePdNOliRpxDNgSdJR6lCWeU+HKVrTu162y8WykOx0Ptl35pBMddDT2k3sZcbiW+KMbZ9EQ9Y6EhOy6E53wpmZx4ul4yS60ky+OIsnbgsJ5seIvxAnfCEkm5wheqU0EPsut7/3+/2/SpIOzoAlSSPYoXS5CEMIAsoYw46wqc+2/sbixOmIt7F1yibCrjCznPsfgQSk16doAx6f00WYFyPWEiM7K5cXS57jxfznCWIBEzqn2M0aprrDLtYWrWJByzmEhKwtWsXJLa/t83VByzkuYCJJL8OAJUkj2CFdzHifixqXBRV9t/U3Rubix4l4nO767sziGg1d9IzqhjnAKZAmID8dMn1sjKeeTcHxAdmtOcSei5H/YpHdrGGgv+5UZ1aS8NQ0W/+yiRAIT03z/INP9vm69S+bGNc9aUDdLTtfko4lnnksSRqQfRfXyGrPZtw9k5m5dT5ZM3KgGLgjIGwJSSWgtT0ktbsbjoP0kynCppBxHRMJguAVH0eDpzvsYm3hKjrC9j5ft5VuIj41i22lW9heuoX41Cy2lzf2+bqtdMtB798Tdh/0sfobk6SRyA6WJGlA+pt2+ETZI+SQQ969BWR1ZbP9ti0k4yHrTuimoCLIrDSYkyD52naWP7SEWBDrXcK9NbZrWHQ3jvYuS3+dpIN1mbZmbyY9p4fHnvgL7dm7uT+rDsZA+9bdcCe0ppszUz+XZI7dub4dqqGzqJ3WSc00jtkIE+D+v9WRmJ9g3f2P93a3snty+zz+vl2xfcckaaSyg3WYkskkDQ0NNDY2DnUpknRE7J1iuP9y8RX3VHPK9rPJSecROyEGlwFzID0PUgnoqegkmBMjcXyCgrGF7EpsJ01qWHQ3jvYuy/4dpZ6w+2W7TNtLtxCfmU1nXpLYSTH4R2BmAPMh6wxgIXA6mWX4PwDB6RB7PcQKIesioCyzb1AWEJuSYGf5NuJTs2gq3nTA4+/tiu07JklHs8bGRhoaGkgmk/2OG7AO0+LFi6murqa2tnaoS5GkIbG3ozW1bRaxIEaYmya7J4fY/43BrdDx/ZDkQyHBqBS5uSEFBT20p9voXtjJ/WV/oGthB5uyX2B7kPlD1d6v+3//aux/zP6OvT1oZGv25kyXJXvTQfcZbC/3GIda//P5T/Y+j73bns2rp3NCB/ePrqM9bzfLR99BY7CRtjt3kUp1k25Nw5+BSZnrnXVvAkYBTUApcDeESyH9fUj/N3R/F3g6hB9DuDuk6/kkndXttBY1s7X0Rdon72b5mDtYMWYJ7fN2k39aIzNO6mB3aTPtr93N8tFLuHf0b1iV/wCdYfJVv7YHe2/7O97LbTuS77GkkaG2tpbq6moWL17c77gB6zB94hOfYOPGjdTV1Q11KZI0JHoXzeCl87PGPTqFqTtOIDYlTnB5jGBBQGoexDpg8sVZBKUBwfzYSx2Q4q0H7cC82o7Hy3Vw9h/rr8vSX1dosLzccz2U7tTe+vc9T2rvOVQ7y7cRn5AFU4HTyHSpXgvBqRB8AGJvgOAi4CkIpgCLyISrEyF4MYBFAezO3I/TyHS1Tt/n/vt0t5gITANmAP8IwWsCWk5MA9BT0RVZB/Ng721/71V/r+1AziGTpP3V1dWxceNGPvGJT/Q77jlYhyk3N5eqqiqampqGuhRJGnJ7u1kllNEU30R2Xi6ddyUJ89KwDtpiIavDrszO+RAbE6O7oJOto1+E4+D+B+rgFFjz0IOUJ8f3WbVuoHq7Ovutgrfpzy+wNX8T7bm7eSDnD4SzIdyVJhFks7u0meSMNpY/lDn5KDEv3uf8oldTx0Bq7e8x9n8evfVn7eb+ijoIoL1nN9wHtEHnne0HnEMFQBrYCkzY06WaAiyH8JQ92ycCncADwEygCcKxwIt7LiKdDTwKnADBRoiNg9TtEO4CUsAYYBuZP91WAr+HsDNkByE7TuyEIiCEVCqkLd1KemGK+//2BxLzE2xa8QITu6cP+PXa/73df6XDg722B7v/YL7HkkaOyspKqqqqDjpuB0uSFJl9u1lbS14kURJn1MZysrbkwlwy3ZFTIHg/5J0VkHtRSN5oYHwAJwFjIZgZo7ugi12l21/qxrzCFLl+vx5kFbzm0m30ZHcROylG+A8QOz1G3kUdHD+nnfzjO3q7LInyRJ/zi/at4+Ued6C2B429tfb3GPs/j731BycHcDKZc6E+CFk1ZDpMr+HAc6jOgNjFmXOpgpY9+zUBC4A1wKYAAojdmyC+OUG8LUHweEC8O07srzGIBQRLY+QU5hE8EcBJAakmMt2q44D3k+lazYfgfZD7VsiaCxwP/HPmfY8tglg3jL8oRqw01qeDuat0+8De44O8t/t28F7utX25FRKdKijpcAVhGIZDXcRI8OijjzJ//nwWcC7FwaihLkeShlQ6TPHwqHup6K6iMFnMutonSGel6FjbDjNCeBzYSebPfCGZDkglL3VAjiPTcZkD3A9sAToDCEImdE9hYuvxrC5+gJNbXstjxX/hhJYFPDb6z6TCHsJEmDlenMx1vHr2HLuJzLyNvf/rVe15vC5e6sA07rdPDdBBJojcu6eOHgjCGPF0nJObX8sTxY/01nEoF+HtDJOsL3qKF3OeByA8Lk1heykshPSKFG3NLZma43vukAIq9tSftU9tU4EWYB6wjkyY6gSeACZD8AKEUyBYD+H8zBS+2POQzoF0855j9ex5viuBAsg+MZfuP3cxZvt4duZsBUKy07lkk0NnOklhqoStb9hE0BEjXdqTuX8nsJ1MlyuE4CnIOhXCJuhOAuv3PNaz+7zfe17bWFuM8LSQ8N4w89ruaW7u+x4f9L1N7zle8NLx2AXkAmv3bJsOtEKsMEZ6bTpzv31/NkIyP2P9vMexIO5FsSX10RLu5GGWsnLlSubNm3fQ/ZwiKEmK3L5TBZ8oe4QgDoV/KWFm0zxWxx4gPTUFrWRCTSeZTsyTEOz5PtjzAToeg+6JAUwLSbyQIFWZIr0iza7s7X2mdjU8sI6gGuJBFj3ju6ALskqgJwnhXOApMlPicsh8iD+ZzIfzx8h8KH/tnn0m7dmnG4J96kjEoOu4AKaHxFsTZD2XTWpSDy8+8NwBU8z2X6Z8/2XSy8IK8nuKyJ2WT+rUHrqKOxldtpvRJXG2pXpo+zuZQJIPBJA4BYK10N3Wt/4gDrQBayA8HVhNpjMFBGsDgpkBwTMB4bwQVqdhYZzU+jTMCAj+FpJzch5dK5OkT09DK+ROyCe9NE3OzFyC3QFzd5xOKuyBIKCMMexiBxvLniU3zCO9Jg0Lsuha2Ul4dgiPQpAOYF4IJ0H3Y3tewxkBpALC5lSmg1mQqT84GXJzAoIgk7ba97y2wfoYVNLnPT7Ye9vd0c/7CZnwXkKmoxaHrBQUbgxpzYfuvExd5EOw7+uY+aGFWZCeDPFkgoLGIvLXelFsSQPnFEFJUuT6W8J9TtMiUokUeZMLiJHIXJS4ESiE4C/AREhMhuDPEDZDVgB5yYD49jjBcTHSO1LkTcgnXZBme2nfqV07xjSRc2oBQWNAbFuCoCBGXnNAYizE7sscm8kQbANGQ7ACgl2QVUqmO7QcKAeqM/skxrxUR3YiIB4EBJuB6QHxp+IU7ComZ0o+28ub+tTxcotk7P2aooexHdWknk2R/lMauqClPbMQxO7NIcGZMRLrskiszYbRmdrikyFnOsS2v1Q/zZBohHAi8MKe59EBPAbhJEhvS5MemyZsSMM4SLenCF8ICZ4KIAwovqeMvN2FjLprDPltRYx7ZFLmwtEPT6ItbKEoLKUsqKCMMQAUhSW0pncx9u6JZLVnM/7uyeTvLmTUneVkt+SQGJsFnRCsiBNOCAm3hMQeiJM9JodgYYysohx4JsisUPhnSDVnwlWYove1DXZwwHt8sPc2axzEdpAJ6SuAFoiFkHMiFIwH/pB5PfKyYkw9M5u8LjJd01woaIbE2MzPXboFSGd+3mKbgAkQezLO7NWnUJ2c5kWxJQ2YUwQj4hRBSTpQT9hNW7C797ysJ8oeYfeZO+l+uBsWQOqvKVKTeghnp2EDMBt4GHiGl6bq7Z0+uHf63t7pc3uneu01B2jf83UlmVXv1pLpnKyG2MJYZkrgMwEcF8IzAekgDdNDYivjkA4JFsVIx1OEj4aZjsZze2oo56Uphj0HedyFQB2ZaWkPAgQk5sUpv388uxZtp/iBMrZN3kx6ayZMhcelKWoN2N0UvvQY24EYxGZArAd6ssmcx/QMme7UNgieCTJhZJ/6gz9nzmFLPJBNqqSH9OwUwaqA2Jw44cNp0meHsAyyxmYTvBgQmxmjZNVoJu48jhLK2MUOCsJi2oKWPl8TwUsTXfa+lwVhUZ+vJZTxWNl9JM9sJ+v+bGKpODvP3Ap/h7ItFbQW7qLr9Z1k/ymH2Q2nsHr2A6TLUpmO4Su9x4f43jIN+AvQvGf/aZDYBj17p6GesOd4kyGrHrpPInPuWXzPcfbsE3QHsDBk8toEmxqz6H42M2fRqYKSwCmCkqRhoL8l3MfcPYGt2ZuouKeK/GQRaxc+RvfmzkzAyCLzIXkuUAgkgVPIfBifTGZKWOd+X/eZIhYEED5LZsGMHcAEoIFMB6crhNUhzApgW5jpfKwCdkKsIkZqTIrsZDZdf+siPLkHniZzXtFCMtPOpvT/uEF8z2lAaQgnQHgSxLbHiO2IE5sSZ+fabWRNzWbn05nl0oPCHsL5IelR0Hl/mmB85hyp3sc4GWIxiG2AIAvCncB4MiFkDYQzQtgWEu6tfzuEXZD1bBbHN8zlmfZVpDfHCIjBspDsdD7Zd2bOoarcWUVT4kUqHp7A1qxNFIWlEJB5j/r52u97ud9+6TBFZ7qDynuqmdR6HI+MWsbku2cA0JT1Ion2LPLvKKQr7KQn0UPutHzaN7fCnDAzlS9JJjD19x7v9xqnD/LesoNMB28ymZCeFZB+MMx0Lk8lMy0wB3gG0rPJBO3JZM4Pq9nzHndB2BMS5MHWGSmyK4vIOjGH2N8CpwpKGhADliTpiNh7XlZBWER5+7je87NyyIHVkFgQp/u2btIladJZezocxwfwpzDzAbiTTAiZTOaD8WQyU+OmkpniNx/CAuBFMl2kVQHBSQE8BOHZIcHtAUFRjHR7iuDxgGBWjFhLQGJCNl1PdsI4CO8OGdVRzs6ZWwm3p6E4gNtCwnOA9jDzeHsfvxL4JYQtmWluvYtkPAHpnjTpeBru6IYUdD7QngkLT5BZXOM+4HV71nQoAe4AziBzXtpy6DkZeGHPeUmrA8KFEPwvBEUB6fZ0n/rjlVmktqTImZvLtuc3cdI+504Rhn3OoSoIiykPxvV+TZMmdphnC+x7vl0P3b3fA4wKxxxwDlcsEZCzOZd0aYruNd2Z87aWkllkor/3eHkmgL7ce8sjIZwVwNKQoCFGOC9NOtjz2v4GeG2QCatbQlLTIVgTEB4PPB9mwtc64Dwy54gByXRIT7gT/gTsCihJ7KaLTnKwgyXplRmwJElHxP4dkP46WtWt03iw/Y+EW+KMbZ/Ei9nPkTe3gPb63YSLQmKzY5kpfrMzU/zCU0LCNWmCU2KZr/NjhA1psnqySe3qofp309mWs4XsO3NIpjroae0mtqe7E9wPi7ZdwIbda1mf9zQld5TRle4kGWtn4u+m82LeeoIdcMq2c1iZt5zu13YSn55F94ZOYrNjhE9DOD2duQjva4EnyXRM8ujbgZlLJiwm6bM4RbiK3nPQOC+AB0Jy5uWSnh7SXd8JL4Zk9WTT091Dzp9zOOWFs3hk9DLSm9N96t/Y+gyNiYbM67inKxUL9oSmfbpQh9KlOqz3db/vAcqCit7v957D1ft+P1TF+NbJPJJ3L12nJqGel3+PD/LeZi3JZldiByV3Zt6/8uZxbAyfIUxmnn/BvEKSf2snd24BbQ27IQwJdgVU/2kam/M3kJ6cInViD8G6gHBRmJmhmAexlXHyzigg/njcDpakATFgSZKGRH8drR66OWnH6QB0xNvYOnkTPU09UBlAJ6RXp1+a4jc2gIYw03nas5BDmExDJcQmx4nHEhTdN4qKjipKKGNH2NSnqxMnTkDApLYZ5LUXUhFW0cIOCAKKw1LGdEwAoCWxg8TULGKdcdIr0zA/IL0tDeOCl84VWk6fZcrj0yD9XGZFuuAve6YAdu5ZnGIexLdDahyE4yHYESN8Nsys7tcdI74SuueQWWxhXYKcE7OJhzF2bdjZpzu1t/6pbbMZzdhIu1KDob/3uymxicTULNKbQ3oqu17+PT7Ie1sQFtEUvEhFWEVb0EJHvI3siXkEWQHBtoCebSnSs0J6urrJLczPhGDSFN9XRmWymviOOKs3PkjygnaCLog/GxI7DrqqQ4K7AmZvPYVCSiIJo5KODcPvN7Ak6Ziwt9vRXwekhDK2l2whEY9DfUhuMo/s27KJtcWhPTzgIrg8TOZiuHfEyM3Kh2WQGJ1ge+mW3mPvXRFv79fMYydIBFmMZzKJINE7tn8dQRwKf1dM4aZi4u3xvo9fGsALAcGaGEHWnovwvhDAjBixDXGYF8ssslESJ5gZI1gdJ7U+gOwY8aeyyH+8gOxVuQSJgPDukIJNxcS7E8SfSFD4ZDFFDxYTJOh9LvvXD/R+v++24aa/93tryYsEcYjVB6/4Hh/svd33/dv7fmXlJCh6sIQT718A9WHmz8nLoGbVQor+WEJidFbv/QspoTg1iqA7IPbnBEF3IT33JAi6oSBVnAlXkjQABixJ0rCz7/TBrPZsxt0zmUUbzyexOUFiaRaJ7dnElsXI211I6ZLRBG0xSu8cTV5zAeMemUSiNYsxd09gd7qZdJg+7Doq7qnmxKYFdKY7qP7dtD6P3+8y5RPDA5dJ70qTXp0mXZXqXcI8d1Y+4WyItQeU3jGa+O4EyXQ7E39/HBOfOY7OWAcnNC2g4u6qw34uw81A3uNRd5Yf0nu7/2UBiigh0Z5F6R2jSbRmUUQpNU2n9Xk902GKXeF24k9n3o94Y5zcWQXEn06wK9w+ol5zSUfG8PwzlyTpmPZy0wf3X8ChOCztM0UsyoUc+lvAoSAsYkzHhD4X4X2s7D6StFOwpIhYKk7zzG2Ej6f7LJPeu8hGa4ogFqP0xdHE6mK0L9xNQVjMybvO6J3GuPfaU+Xt4wB6pwEO1+l/r8ZA3uNDXaRj3/eLAFJhmtnNp/QuRZ8mTSJI9Hk9IaQgXUznxHby6gqY3DST9c88RXBGSM66/BH1mks6MrwOVkS8DpYkHZvSYYqHR91LRXdV7zLl5Z1jeTFvPRCSnc4lmz2LbGR1QxAyoWMq27I3c8rO17Gh8Gmasl5kwc5zX1qcQkfMvu/flNaZBEFAOkyzvvAp3xdJfXgdLEmSjoBD7XLt7U7FiVMQFjOmfTwwMrtTR5P9u14AsSDm+yLpVTNgSZJ0GPpdpryfZdD3XbJ8//GolkzXwO2/tPy+fF8kvRr+SUaSJEmSImLAkiRJkqSIGLAkSZIkKSIGLEmSJEmKiAFLkiRJkiJiwJIkSZKkiBiwJEmSJCkiBixJkiRJiogBS5IkSZIiYsCSJEmSpIgYsCRJkiQpIgYsSZIkSYrIiAxYN954I3PmzKG4uJji4mIWLVrEH/7wh97xZDLJVVddxejRoyksLOTSSy+lsbGxd3z9+vUEQTAUpUuSJEk6io3IgFVVVcW///u/s3LlSv72t79xzjnn8KY3vYknnngCgKuvvpo777yT22+/nRUrVrBp0ybe+ta3DnHVkiRJko52iaEuYDBcfPHFfW5//etf58Ybb+Shhx6iqqqKH/3oR9x6662cc845ANx8883MmjWLhx56iNNOO20oSpYkSZI0AozIDta+UqkUv/jFL2hra2PRokWsXLmS7u5uzjvvvN59Zs6cycSJE3nwwQeHsFJJkiRJR7sR2cECqK+vZ9GiRSSTSQoLC/nNb37D7NmzWbVqFdnZ2ZSWlvbZv7Kyki1btgAwefJkwjAcgqolSZIkHc1GbMCaMWMGq1atYteuXfy///f/uOyyy1ixYkVkx29tbaW5ubn39r6LZEiSJEkamRobG2loaOi9XVpaSmFhYe/tERuwsrOzmT59OgDz58/nkUce4b/+679429veRldXF83NzX26WI2NjYwdO/aQj/+pT32Km266KeqyJUmSJA1jtbW1fW5/8IMf5MYbb+y9PeLPwdornU7T2dnJ/PnzycrKYunSpb1ja9euZcOGDSxatOiQj/etb32LjRs39v6rq6sbjLIlSZIkDSN1dXV9csC3vvWtPuMjsoP12c9+lgsvvJCJEyeye/dubr31VpYvX84f//hHSkpKeN/73sc111xDWVkZxcXFfOQjH2HRokUDWkGwsLCwTyuwqalpMJ6KJEmSpGGksrKSqqqqg46PyIDV1NTEe97zHjZv3kxJSQlz5szhj3/8I69//esBuP7664nFYlx66aV0dnZy/vnn8/3vf3+Iq5YkSZJ0tAtCl8uLxKOPPsr8+fNZwLkUB6OGuhxJkiRJEWoJd/IwS1m5ciXz5s076H7HzDlYkiRJkjTYDFiSJEmSFBEDliRJkiRFxIAlSZIkSRExYEmSJElSRAxYkiRJkhQRA5YkSZIkRcSAJUmSJEkRMWBJkiRJUkQMWJIkSZIUEQOWJEmSJEXEgCVJkiRJETFgSZIkSVJEDFiSJEmSFBEDliRJkiRFxIAlSZIkSRExYEmSJElSRAxYkiRJkhQRA5YkSZIkRcSAJUmSJEkRMWBJkiRJUkQMWJIkSZIUEQOWJEmSJEXEgCVJkiRJETFgSZIkSVJEDFiSJEmSFBEDliRJkiRFxIAlSZIkSRExYEmSJElSRAxYhymZTNLQ0EBjY+NQlyJJkiRpkDU2NtLQ0EAymex33IB1mBYvXkx1dTW1tbVDXYokSZKkQVZbW0t1dTWLFy/udzwIwzA8wjWNKMlkkm3btlFfX09tbS0LOJfiYNRQlyVJkiQpQi3hTh5mKXV1ddTU1FBeXk5ubu4B+yWGoLYRJTc3l6qqKpqamoa6FEmSJEmDrLKykqqqqoOOO0VQkiRJkiJiwJIkSZKkiBiwJEmSJCkiBixJkiRJiogBS5IkSZIiYsCSJEmSpIgYsCRJkiQpIgYsSZIkSYqIAUuSJEmSImLAkiRJkqSIGLAkSZIkKSIGLEmSJEmKiAFLkiRJkiJiwJIkSZKkiBiwJEmSJCkiBixJkiRJiogBS5IkSZIiYsCSJEmSpIgYsCRJkiQpIgYsSZIkSYqIAUuSJEmSImLAkiRJkqSIGLAkSZIkKSIGLEmSJEmKiAFLkiRJkiJiwJIkSZKkiBiwJEmSJCkiBixJkiRJiogBS5IkSZIiYsCSJEmSpIgYsCRJkiQpIgYsSZIkSYqIAeswJZNJGhoaaGxsHOpSJEmSJA2yxsZGGhoaSCaT/Y4bsA7T4sWLqa6upra2dqhLkSRJkjTIamtrqa6uZvHixf2OB2EYhke4phElmUyybds26uvrqa2tZQHnUhyMGuqyJEmSJEWoJdzJwyylrq6OmpoaysvLyc3NPWC/xBDUNqLk5uZSVVVFU1PTUJciSZIkaZBVVlZSVVV10HGnCEqSJElSRAxYkiRJkhQRA5YkSZIkRcSAJUmSJEkRMWBJkiRJUkQMWJIkSZIUEQOWJEmSJEXEgCVJkiRJETFgSZIkSVJEDFiSJEmSFBEDliRJkiRFxIAlSZIkSRExYEmSJElSRAxYkiRJkhQRA5YkSZIkRcSAJUmSJEkRMWBJkiRJUkQMWJIkSZIUEQOWJEmSJEXEgCVJkiRJETFgSZIkSVJEDFiSJEmSFBEDliRJkiRFxIAlSZIkSRExYEmSJElSRAxYkiRJkhQRA5YkSZIkRcSAJUmSJEkRMWBJkiRJUkQMWJIkSZIUEQOWJEmSJEXEgCVJkiRJETFgHaZkMklDQwONjY1DXYokSZKkQdbY2EhDQwPJZLLfcQPWYVq8eDHV1dXU1tYOdSmSJEmSBlltbS3V1dUsXry43/EgDMPwCNc0oiSTSbZt20Z9fT21tbUs4FyKg1FDXZYkSZKkCLWEO3mYpdTV1VFTU0N5eTm5ubkH7JcYgtpGlNzcXKqqqmhqahrqUiRJkiQNssrKSqqqqg467hRBSZIkSYqIAUuSJEmSImLAkiRJkqSIGLAkSZIkKSIGLEmSJEmKiAFLkiRJkiJiwJIkSZKkiBiwJEmSJCkiBixJkiRJiogBS5IkSZIiYsCSJEmSpIgYsCRJkiQpIgYsSZIkSYqIAUuSJEmSImLAkiRJkqSIGLAkSZIkKSIGLEmSJEmKiAFLkiRJkiJiwJIkSZKkiBiwJEmSJCkiBixJkiRJiogBS5IkSZIiYsCSJEmSpIgYsCRJkiQpIgYsSZIkSYqIAUuSJEmSImLAkiRJkqSIGLAkSZIkKSIGLEmSJEmKiAFLkiRJkiJiwJIkSZKkiBiwJEmSJCkiBqzDlEwmaWhooLGxcahLkSRJkjTIGhsbaWhoIJlM9jtuwDpMixcvprq6mtra2qEuRZIkSdIgq62tpbq6msWLF/c7HoRhGB7hmkaUZDLJtm3bqK+vp7a2lgWcS3EwaqjLkiRJkhShlnAnD7OUuro6ampqKC8vJzc394D9EkNQ24iSm5tLVVUVTU1NQ12KJEmSpEFWWVlJVVXVQcedIihJkiRJETFgSZIkSVJEDFiSJEmSFBEDliRJkiRFxIAlSZIkSRExYEmSJElSRAxYkiRJkhQRA5YkSZIkRcSAJUmSJEkRMWBJkiRJUkQMWJIkSZIUEQOWJEmSJEXEgCVJkiRJETFgSZIkSVJEDFiSJEmSFBEDliRJkiRFxIAlSZIkSRExYEmSJElSRAxYkiRJkhQRA5YkSZIkRcSAJUmSJEkRMWBJkiRJUkQMWJIkSZIUEQOWJEmSJEXEgCVJkiRJETFgSZIkSVJEDFiSJEmSFJERH7D+/d//nSAI+PjHP9677eyzzyYIgj7/PvjBD/aOr1+/niAIhqBaSZIkSUezxFAXMJgeeeQRfvCDHzBnzpwDxt7//vfz1a9+tfd2fn7+kSxNkiRJ0gg0YjtYra2tvOtd7+J//ud/GDVq1AHj+fn5jB07tvdfcXHxEFQpSZIkaSQZsQHrqquu4qKLLuK8887rd/znP/855eXlnHjiiXz2s5+lvb39CFcoSZIkaaQZkVMEf/GLX/Doo4/yyCOP9Dv+zne+k0mTJjF+/HjWrFnDZz7zGdauXcuvf/1rACZPnkwYhkeyZEmSJEkjwIgLWBs3buRjH/sY99xzD7m5uf3u84EPfKD3+5qaGsaNG8e5557LunXrmDZt2iE9TmtrK83Nzb23GxsbD6tuSZIkScNfY2MjDQ0NvbdLS0spLCzsvR2EI6xV89vf/pa3vOUtxOPx3m2pVIogCIjFYnR2dvYZA2hra6OwsJC77rqL888//5Ae58orr+Smm246YPsCzqU4OPCcL0mSJElHr5ZwJw+z9IDtH/zgB7nxxht7b4+4Dta5555LfX19n21XXHEFM2fO5DOf+cwB4Qpg1apVAIwbN+6QH+db3/oWn/vc53pv19fXU1tb++qKliRJknRUqKuro6ampvd2aWlpn/ERF7CKioo48cQT+2wrKChg9OjRnHjiiaxbt45bb72V2tpaRo8ezZo1a7j66qs588wz+13O/WAKCwv7tAKbmpoiew6SJEmShqfKykqqqqoOOj7iAtYryc7O5k9/+hM33HADbW1tVFdXc+mll/L5z39+qEuTJEmSdJQ7JgLW8uXLe7+vrq5mxYoVQ1eMJEmSpBFrxF4HS5IkSZKONAOWJEmSJEXEgCVJkiRJETFgSZIkSVJEDFiSJEmSFBEDliRJkiRFxIAlSZIkSRExYEmSJElSRAxYkiRJkhQRA5YkSZIkRcSAJUmSJEkRMWBJkiRJUkQMWJIkSZIUEQOWJEmSJEXEgCVJkiRJETFgSZIkSVJEDFiSJEmSFBEDliRJkiRFxIAlSZIkSRExYEmSJElSRAxYkiRJkhQRA5YkSZIkRcSAJUmSJEkRMWBJkiRJUkQMWJIkSZIUEQOWJEmSJEXEgCVJkiRJETFgSZIkSVJEDFiSJEmSFBED1mFKJpM0NDTQ2Ng41KVIkiRJGmSNjY00NDSQTCb7HTdgHabFixdTXV1NbW3tUJciSZIkaZDV1tZSXV3N4sWL+x0PwjAMj3BNI0oymWTbtm3U19dTW1vLAs6lOBg11GVJkiRJilBLuJOHWUpdXR01NTWUl5eTm5t7wH6JIahtRMnNzaWqqoqmpqahLkWSJEnSIKusrKSqquqg404RlCRJkqSIGLAkSZIkKSIGLEmSJEmKiAFLkiRJkiJiwJIkSZKkiBiwJEmSJCkiBixJkiRJiogBS5IkSZIiYsCSJEmSpIgYsCRJkiQpIgYsSZIkSYqIAUuSJEmSImLAkiRJkqSIGLAkSZIkKSIGLEmSJEmKiAFLkiRJkiJiwJIkSZKkiBiwJEmSJCkiBixJkiRJiogBS5IkSZIiYsCSJEmSpIgYsCRJkiQpIgYsSZIkSYqIAUuSJEmSIpIY6gJGio6ODgA++rMrmDVr1hBXI0mSJClKTz75JO9+99Lez/0HY8CKyPr16wF497vfPbSFSJIkSRo069ev5/TTTz/oeBCGYXgE6xmxtm3bxh//+EcmT55MXl7eoD1OY2MjtbW11NXVUVlZOWiPczTWM5xqsR7rGSm1WI/1jJRarMd6Rkot1jN09XR0dLB+/XrOP/98ysvLD7qfHayIlJeX8653vWvQH6ehoQGAmpoaqqqqBv3xXslwqmc41QLW80qs5+ioBaznlVjP0VELWM8rsZ6joxawnlcymPW8XOdqLxe5OMqUl5fzta997WVT85E0nOoZTrUMR8Pt9SkuLuass86iuLh4qEsBhlc9w6mW4Wi4vT7Wc3TUMhwNt9fHeo6OWoaj4fYZY6jrcYqgNEgaGhqorq5m48aNw+KvOdKr5c+yRgp/ljVS+LM8vNnBkgZJaWkpH/zgByktLR3qUqTD4l9uNVL4s6yRws8Yw5sdLEmSJEmKiB0sSZIkSYqIAUuSJEmSImLAkiRJkqSIGLAkSZIkKSIGLAn45je/yamnnkpRUREVFRW8+c1vZu3atX32OfvsswmCoM+/D37wg332WbJkCccffzwzZszgd7/7Xe/27du3c8EFFzB+/HhycnKorq7mwx/+MC0tLX3uv3z5cubNm0dOTg7Tp0/nlltu6TN++eWX8+UvfznS5y79+c9/5uKLL2b8+PEEQcBvf/vbPuOXX375AT/7F1xwQZ99giBg/fr1R65o6SC+973vMXnyZHJzc1m4cCEPP/xw71gymeSqq65i9OjRFBYWcumll9LY2Ng7vn79eoIgGIqyNYIdymeMvcIw5MILL+z3d/HBPmPsa/v27VRVVREEAc3Nzb3bly9ffsDv8SAI2LJlS+8+fsaIjgFLAlasWMFVV13FQw89xD333EN3dzdveMMbaGtr67Pf+9//fjZv3tz77z//8z97xzo7O7nqqqv4/ve/z3e/+12uvPJKurq6AIjFYrzpTW9iyZIlPP3009xyyy386U9/6hPQnn/+eS666CJe97rXsWrVKj7+8Y/zz//8z/zxj388Mi+CjlltbW2cdNJJfO973zvoPhdccEGfn/3//d//PYIVSofmtttu45prruFLX/oSjz76KCeddBLnn38+TU1NAFx99dXceeed3H777axYsYJNmzbx1re+dYir1kh3qJ8xAG644YZ+Q/7LfcbY1/ve9z7mzJlz0FrWrl3b53d5RUXF4T059Ssx1AVIw8Fdd93V5/Ytt9xCRUUFK1eu5Mwzz+zdnp+fz9ixY/s9RmdnJ/F4nLlz5wKQSCTo7OwkOzubUaNGceWVV/buO2nSJD70oQ/xrW99q3fbTTfdxJQpU1i8eDEAs2bN4r777uP666/n/PPPj+qpSge48MILufDCC192n5ycnIP+7EvDxXXXXcf73/9+rrjiCiDze/X3v/89P/7xj7nyyiv50Y9+xK233so555wDwM0338ysWbN46KGHOO2004aydI1gh/oZY9WqVSxevJi//e1vjBs3rs99Xu4zxl433ngjzc3NfPGLX+QPf/hDv7VUVFR47awjwA6W1I9du3YBUFZW1mf7z3/+c8rLyznxxBP57Gc/S3t7e+9YcXExV1xxBePGjWP8+PFceeWVFBUV9Xv8TZs28etf/5qzzjqrd9uDDz7Ieeed12e/888/nwcffDCqpyW9asuXL6eiooIZM2Zw5ZVXsn379qEuSeqjq6uLlStX9vk9GovFOO+883jwwQdZuXIl3d3dfcZnzpzJxIkT/T2rI6q/zxjt7e28853v5Hvf+16/f8x6pc8Yf//73/nqV7/KT3/6U2Kxg3+8nzt3LuPGjeP1r389999/f4TPSvuygyXtJ51O8/GPf5zTTz+dE088sXf7O9/5TiZNmsT48eNZs2YNn/nMZ1i7di2//vWve/f50pe+xMc//nFisVi/4eod73gHd9xxBx0dHVx88cX88Ic/7B3bsmULlZWVffavrKykpaWFjo4O8vLyDjgnSzoSLrjgAt761rcyZcoU1q1bx7/9279x4YUX8uCDDxKPx4HMeQPSUNq2bRupVKrf36NPPfUUW7ZsITs7+4C/3ldWVvaehzJ58mR/ljWoDvYZ4+qrr+Y1r3kNb3rTmw5634N9xujs7OQd73gH3/rWt5g4cSLPPffcAfcdN24cN910E6eccgqdnZ388Ic/5Oyzz+avf/0r8+bNA/AzRoQMWNJ+rrrqKh5//HHuu+++Pts/8IEP9H5fU1PDuHHjOPfcc1m3bh3Tpk3rHSspKTnosa+//nq+9KUv8fTTT/PZz36Wa665hu9///vRPwkpQm9/+9t7v6+pqWHOnDlMmzaN5cuXc+655w5hZZJ0dOnvM8aSJUu49957eeyxx17x/v19xvjsZz/LrFmzePe7333Q+82YMYMZM2b03n7Na17DunXruP766/m///f/DvBZ6JU4RVDax4c//GF+97vfsWzZMqqqql5234ULFwLw7LPPHvLxx44dy8yZM7nkkkv4wQ9+wI033sjmzZt7x/ZdzQqgsbGR4uJi8vLyBvhMpMEzdepUysvLB/SzLw228vJy4vF4v79Hx44dy9ixY+nq6uqzstq+49JgO9hnjHvvvZd169ZRWlpKIpEgkcj0Py699FLOPvvsVzzuvffey+233957371/+CovL+dLX/rSQe+3YMECf48PEjtYEpnpTR/5yEf4zW9+w/Lly5kyZcor3mfVqlUAB5yIeqjS6TSQae0DLFq0iLq6uj773HPPPSxatOhVHV8aLA0NDWzfvv1V/+xLgyE7O5v58+ezdOlS3vzmNwOZ37NLly7lwx/+MPPnzycrK4ulS5dy6aWXApkV1TZs2ODvWQ2qV/qM8a//+q/88z//c59tNTU1XH/99Vx88cWvePxf/epXdHR09N5+5JFHeO9738tf/vKXPjNs9rdq1Sp/jw8SA5ZEpmV/6623cscdd1BUVNQ7H7+kpIS8vDzWrVvHrbfeSm1tLaNHj2bNmjVcffXVnHnmmS+7HOpedXV1NDY2cuqpp1JYWMgTTzzBpz71KU4//XQmT54MwAc/+EG++93v8ulPf5r3vve93Hvvvfzyl7/k97///WA+dYnW1tY+f8V8/vnnWbVqFWVlZZSVlfGVr3yFSy+9lLFjx7Ju3To+/elPM336dFe31LBzzTXXcNlll3HKKaewYMECbrjhBtra2rjiiisoKSnhfe97H9dccw1lZWUUFxfzkY98hEWLFrmCoAbVK33G2Nth3d/EiRMP6Q+++4eobdu2AZnViPeec3jDDTcwZcoUTjjhBJLJJD/84Q+59957ufvuuw/z2alfoaQQ6PffzTffHIZhGG7YsCE888wzw7KysjAnJyecPn16+KlPfSrctWvXIR3/3nvvDRctWhSWlJSEubm54XHHHRd+5jOfCXfu3Nlnv2XLloVz584Ns7Ozw6lTp/Y+vjSYli1b1u/P/2WXXRa2t7eHb3jDG8IxY8aEWVlZ4aRJk8L3v//94ZYtW4a6bKlf3/nOd8KJEyeG2dnZ4YIFC8KHHnqod6yjoyP80Ic+FI4aNSrMz88P3/KWt4SbN28ewmp1LHilzxgHu89vfvObV/V4e3+n7/sZ4z/+4z/CadOmhbm5uWFZWVl49tlnh/fee++rOr5eWRCGLpcjSZIkSVFwkQtJkiRJiogBS5IkSZIiYsCSJEmSpIgYsCRJkiQpIgYsSZIkSYqIAUuSJEmSImLAkiRJkqSIGLAkSZIkKSIGLEnSsHDLLbdQWlo66I9z9tln8/GPf3zQH2cwLV++nCAIaG5ujvzY69evJwgCVq1aFfmxJelYYMCSJB2WLVu28LGPfYzp06eTm5tLZWUlp59+OjfeeCPt7e1DVtfBQsivf/1rrr322ld93OEQ0F7zmtewefNmSkpKgCMXTiVJrywx1AVIko5ezz33HKeffjqlpaV84xvfoKamhpycHOrr6/nv//5vJkyYwCWXXDLUZfZRVlY21CUctuzsbMaOHTvUZUiS+mEHS5L0qn3oQx8ikUjwt7/9jX/8x39k1qxZTJ06lTe96U38/ve/5+KLL+7d97rrrqOmpoaCggKqq6v50Ic+RGtr60GPvW7dOt70pjdRWVlJYWEhp556Kn/605/67NPZ2clnPvMZqqurycnJYfr06fzoRz9i/fr1vO51rwNg1KhRBEHA5ZdfDhzYgTrYMV6tX/3qV5xwwgnk5OQwefJkFi9e3Gd88uTJfOMb3+C9730vRUVFTJw4kf/+7//us88DDzzA3Llzyc3N5ZRTTuG3v/1tn2l7+3bnli9fzhVXXMGuXbsIgoAgCPjyl78MQBAE/Pa3v+1z7NLSUm655Zbe2w8//DAnn3xy72M99thjBzynxx9/nAsvvJDCwkIqKyv5p3/6J7Zt2/aqXyNJGskMWJKkV2X79u3cfffdXHXVVRQUFPS7TxAEvd/HYjG+/e1v88QTT/CTn/yEe++9l09/+tMHPX5rayu1tbUsXbqUxx57jAsuuICLL76YDRs29O7znve8h//93//l29/+Nk8++SQ/+MEPKCwspLq6ml/96lcArF27ls2bN/Nf//Vf/T7OwY7xaqxcuZJ//Md/5O1vfzv19fV8+ctf5gtf+EKfQAOwePHi3jDzoQ99iCuvvJK1a9cC0NLSwsUXX0xNTQ2PPvoo1157LZ/5zGcO+pivec1ruOGGGyguLmbz5s1s3ryZT37yk4dUb2trK2984xuZPXs2K1eu5Mtf/vIB921ubuacc87h5JNP5m9/+xt33XUXjY2N/OM//uPAXhxJOlaEkiS9Cg899FAIhL/+9a/7bB89enRYUFAQFhQUhJ/+9KcPev/bb789HD16dO/tm2++OSwpKXnZxzzhhBPC73znO2EYhuHatWtDILznnnv63XfZsmUhEO7cubPP9rPOOiv82Mc+dkjH6M++99/fO9/5zvD1r399n22f+tSnwtmzZ/fenjRpUvjud7+793Y6nQ4rKirCG2+8MQzDMLzxxhvD0aNHhx0dHb37/M///E8IhI899li/z+1grx0Q/uY3v+mzraSkJLz55pvDMAzDH/zgBwc81o033tjnsa699trwDW94Q59jbNy4MQTCtWvX9vs6SNKxzA6WJClSDz/8MKtWreKEE06gs7Ozd/uf/vQnzj33XCZMmEBRURH/9E//xPbt2w+6EEZrayuf/OQnmTVrFqWlpRQWFvLkk0/2drBWrVpFPB7nrLPOetW1RnGMfT355JOcfvrpfbadfvrpPPPMM6RSqd5tc+bM6f0+CALGjh1LU1MTkOm4zZkzh9zc3N59FixYEEl9/dW7/2MtWrSozz6rV69m2bJlFBYW9v6bOXMmkJnGKUnqy0UuJEmvyvTp0wmCoHdq215Tp04FIC8vr3fb+vXreeMb38iVV17J17/+dcrKyrjvvvt43/veR1dXF/n5+Qcc/5Of/CT33HMP/+f//B+mT59OXl4e/9//9//R1dV1wPFfrSiO8WpkZWX1uR0EAel0OvLHCYKAMAz7bOvu7h7QMVpbW7n44ov5j//4jwPGxo0bd1j1SdJIZAdLkvSqjB49mte//vV897vfpa2t7WX3XblyJel0msWLF3Paaadx/PHHs2nTppe9z/3338/ll1/OW97yFmpqahg7dizr16/vHa+pqSGdTrNixYp+75+dnQ3Qp3O0v1c6xkDNmjWL+++/v8+2+++/n+OPP554PH5Ix5gxYwb19fV9un+PPPLIy94nOzu73+c5ZswYNm/e3Hv7mWee6dMxnDVrFmvWrCGZTPZue+ihh/ocY968eTzxxBNMnjyZ6dOn9/l3sHPvJOlYZsCS/v/27icUmjiO4/h3Dju1cVAuJrFot+RgFbmKA2cHyl42irY4kIMbZXNacbDJSVl/ipCDtvaw5SIHyu7Fri32RKuVEzmQ73N6tta/R09TPD3vV81h/v3mz+3TzHwGwF9bXFyU5+dnaWlpkc3NTUmlUnJ+fi5ra2uSTqcLocLtdsvT05MsLCzI5eWlrK6uytLS0qdjezwe2d3dlUQiIclkUnw+X9FTnpqaGvH7/TIwMCB7e3uSzWbl4OBAtra2RETE5XKJYRiyv78v+Xz+3cbCP43xkXw+L4lEomi6ubmR8fFxicfjEgwGJZPJyMrKioTD4S+XTohI4TqHhoYklUpJLBaT2dlZESkuDXl9Hff39xKPx+X29rYQojo6OiQcDsvp6amcnJxIIBAoenrm8/nEMAwZHByUs7MziUajhWP9Njw8LHd3d9LX1yfHx8dycXEhsVhM+vv7Pw2vAPDf+u6PwAAA/7br62sdGRnR2tpadTgcWlpaqq2trRoKhfTh4aGw3dzcnFqWpU6nU7u6ujQSiXxa1JDNZrW9vV2dTqdWVVVpOBx+UzDx+PioY2NjalmWmqapbrdbl5eXC+unp6e1oqJCDcNQv9+vqm9LKv40xmttbW0qIm+mYDCoqqrb29va0NCgDodDq6urNRQKFe3vcrl0fn6+aJnX69WpqanC/OHhoTY2Nqppmtrc3KwbGxsqIppOp1X1/QKPQCCg5eXlKiKFsa6urrSzs1NLSkrU4/FoNBotKrlQVT06OlKv16umaWpTU5Pu7OwUlVyoqmYyGe3u7taysjJ1Op1aX1+vo6Oj+vLy8uF9AoD/laH66uVsAADwo6yvrxf+dfVd340BAL6GkgsAAH6YSCQidXV1UllZKclkUiYmJqS3t5dwBQD/AAIWAAA/TC6Xk8nJScnlcmJZlvT09MjMzMx3nxYA4At4RRAAAAAAbEKLIAAAAADYhIAFAAAAADYhYAEAAACATQhYAAAAAGATAhYAAAAA2ISABQAAAAA2IWABAAAAgE0IWAAAAABgEwIWAAAAANjkF4uk6zH8t20eAAAAAElFTkSuQmCC",
      "text/plain": [
       "<Figure size 1000x900 with 1 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAj4AAAGdCAYAAAASUnlxAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjYuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/P9b71AAAACXBIWXMAAA9hAAAPYQGoP6dpAABU5ElEQVR4nO3de1yUdfo//tc9gEAqjCgIJBpqhkRpoiJlmcIaYmZ920zC7eRKmdZmbh9lazPbWvXjbtrR1jRbf5i19umEu9KimJqSKEqGICriCZlMRiDkzLx/f0wzAc6Zuef4ej4e89iduU/vm7mRq/t+X9clCSEEiIiIiLyAwtkDICIiInIUBj5ERETkNRj4EBERkddg4ENEREReg4EPEREReQ0GPkREROQ1GPgQERGR12DgQ0RERF7D19kDcAaNRoMLFy6gd+/ekCTJ2cMhIiIiCwgh8PPPPyMyMhIKhW33brwy8Llw4QKioqKcPQwiIiKywblz5zBgwACbtvXKwKd3794AgGvnfgiF/zX6zz94ZAzGDg5x1rCIiIjIhLq6OkRFRen/jtvCKwMf3eMthf81+sBHGeiHwF690LNXb/go+PiLiIjIVXVnmgonN/+iprEV6ev2Y/yKPOQUVzl7OERERCQDBj5dqGqbMDfrEIMfIiIiD8TApwvxy/8uzS5Bu0aYXJeIiIjcCwMfAwSAqtomFFSonT0UIiIisiMGPiZc/LnJ2UMgIiIiO2LgY0JY7wBnD4GIiIjsyKsDH1PJcCE9/RA/qI/DxkJERETy8+rABzAe/KivtGLCyp3M7iIiIvIgXh34vP7gCIQHG3+cxdR2IiIiz+LVgc9vYsOx6/mJCOnZw+ByprYTERF5Fq8OfACg8MxlqK+0GF3O1HYiIiLP4fWBj6Up60xtJyIicn9eH/hYmrLO1HYiIiL35/WBz9joEEQEBxjN7pIARAQHYGx0iCOHRURERDLw+sDHRyFhybRYAIZT2wWA1LhwFFSoOcGZiIjIzUlCCK/7a15XV4fg4GDU1tYiKCgIAJBTXIWl2SWoqv11Lo9CAjrGOhHBAVgyLRYpcRGOHjIREZHXM/T321oMfDr84No1AgUVauSWqPDB3tNXbae7I7Rm1igGP0RERA5mj8DH6x91deSjkDA2OgTbilUGl7OuDxERkXtj4NNFQYW60+OurljXh4iIyH0x8OmCdX2IiIg8FwOfLljXh4iIyHMx8OmCdX2IiIg8FwOfLnR1fYxNXRYAlkyLhY/CWGhEREREroqBDxEREXkNBj5dtGsElmaXGF0ugensRERE7krWwEetViM9PR1BQUFQKpWYPXs26uvrTW7zxBNPYMiQIQgMDERoaCimT5+OY8eO6Zd///33SEtLQ1RUFAIDAzF8+HC88cYbdhsz09mJiIg8l6yBT3p6Oo4ePYrc3Fxs3boVu3fvRkZGhslt4uPjsWHDBpSWluLrr7+GEAKTJ09Ge3s7AKCwsBBhYWHIysrC0aNH8cILLyAzMxNvv/22XcbMdHYiIiLPJVvLitLSUsTGxuLAgQMYPXo0ACAnJwepqak4f/48IiMjLdrPkSNHMGLECJw8eRJDhgwxuM68efNQWlqKvLw8i/ZpquR1fnk10t7/zuw+Ns8Zh8QhfS06HhEREXWfS7esyM/Ph1Kp1Ac9AJCcnAyFQoH9+/dbtI8rV65gw4YNiI6ORlRUlNH1amtrERJiPL28ubkZdXV1nV7GMJ2diIjIc8kW+KhUKoSFhXX6zNfXFyEhIVCpDPfC0nn33XfRq1cv9OrVC9u2bUNubi569OhhcN19+/bhk08+MfkIbdmyZQgODta/TAVRunR2AAaDHwFg5hjj2xMREZHrsjrwWbx4MSRJMvnqOBnZFunp6Th8+DB27dqFYcOGYcaMGWhqunpOTXFxMaZPn44lS5Zg8uTJRveXmZmJ2tpa/evcuXMmj58SF4E1s0YhPNhwdeZV209g/Io85BRXWXdiRERE5FRWz/H56aefUF1dbXKdwYMHIysrCwsXLsTly5f1n7e1tSEgIABbtmzBfffdZ9HxWlpa0KdPH6xbtw5paWn6z0tKSjBx4kT8/ve/x2uvvWbNKVj8jLBdI/B23kms2n78qmW6u0FrZo1CSlyEVccnIiIi69ljjo+vtRuEhoYiNDTU7HqJiYmoqalBYWEh4uPjAQB5eXnQaDRISEiw+HhCCAgh0NzcrP/s6NGjmDRpEh555BGrgx5rfXzgrOFx4deaPr+JDWclZyIiIjcg2xyf4cOHIyUlBXPmzEFBQQH27t2L+fPnY+bMmfqMrsrKSsTExKCgoAAAcOrUKSxbtgyFhYU4e/Ys9u3bhwceeACBgYFITU0FoH28NXHiREyePBnPPfccVCoVVCoVfvrpJ7ufA2v6EBEReRZZ6/hs2rQJMTExSEpKQmpqKsaPH4+1a9fql7e2tqKsrAwNDQ0AgICAAOzZswepqakYOnQoHnzwQfTu3Rv79u3TT5T+9NNP8dNPPyErKwsRERH615gxY+w+ftb0ISIi8iyy1fFxZZY+I2RNHyIiItfh0nV8PIG5mj4AoAz0g0YI9u4iIiJyAwx8TDBX0wcAahpbkb5uP9PbiYiI3AADHzPM1fTRUdU2YW7WIQY/RERELoyBjwVS4iLw7aJJ2DQ7AcpAP4Pr6B50Lc0u4WMvIiIiF8XAx0I+CgkKhYSaxlaj6zC9nYiIyLUx8LEC09uJiIjcGwMfK4T1Nj3Px9r1iIiIyLEY+FiB6e1ERETujYGPFZjeTkRE5N4Y+FiJ6e1ERETui4GPDZjeTkRE5J4Y+NiI6e1ERETuh4FPNzC9nYiIyL0w8OkGprcTERG5FwY+3WAuvV0CEBEcgLHRIY4cFhERERnBwKcbTKW3694vmRYLH4Wpyj9ERETkKAx8uiklLgIZd0RD6hLbSBKQcUc0UuIinDMwIiIiugoDn27KKa7C2t0V6JqxrhHA2t0VrONDRETkQhj4dEO7RmBpdglMVelhHR8iIiLXwcCnGwoq1KiqNZ6qzjo+REREroWBTzewjg8REZF7YeDTDazjQ0RE5F4Y+HSDuTo+ABDS0w/xg/o4bExERERkHAOfbjBVx0dHfaUVE1buZHYXERGRC2Dg000pcRFYM2sUwoONP85S1TZhbtYhBj9EREROxsDHDlLiIrDr+YkI6dnD4HJdMjtT24mIiJyLgY+dFJ65DPWVFqPLmdpORETkfAx87ISp7URERK6PgY+dMLWdiIjI9THwsRNzqe0SgIjgAIyNDnHksIiIiKgDBj52Yiq1Xfd+ybRY+ChMVf0hIiIiOTHwsaOUuAhk3BENqUtsI0lAxh3RSImLcM7AiIiICAADH7vKKa7C2t0V6JqxrhHA2t0VrONDRETkZAx87KRdI7A0uwSmqvSwjg8REZFzMfCxk4IKNapqjaeqs44PERGR88ka+KjVaqSnpyMoKAhKpRKzZ89GfX29yW2eeOIJDBkyBIGBgQgNDcX06dNx7Ngxg+tWV1djwIABkCQJNTU1MpyB5VjHh4iIyPXJGvikp6fj6NGjyM3NxdatW7F7925kZGSY3CY+Ph4bNmxAaWkpvv76awghMHnyZLS3t1+17uzZs3HzzTfLNXyrsI4PERGR65Mt8CktLUVOTg7WrVuHhIQEjB8/Hm+99RY+/vhjXLhwweh2GRkZuOOOO3Dddddh1KhRePXVV3Hu3DmcPn2603pr1qxBTU0N/vjHP8p1ClYxV8cHAJSBftAIwXk+RERETiJb4JOfnw+lUonRo0frP0tOToZCocD+/fst2seVK1ewYcMGREdHIyoqSv95SUkJXnnlFWzcuBEKhWtMUzJVx0enprEV6ev2Y/yKPGZ4EREROYFsUYNKpUJYWFinz3x9fRESEgKVSmVy23fffRe9evVCr169sG3bNuTm5qJHD23n8+bmZqSlpWHlypUYOHCgRWNpbm5GXV1dp5ccUuIisGbWKIQHm36cpaptwtysQwx+iIiIHMzqwGfx4sWQJMnky9hkZEulp6fj8OHD2LVrF4YNG4YZM2agqUk7KTgzMxPDhw/HrFmzLN7fsmXLEBwcrH91vHtkbylxEfh20SRsmp0AZaCfwXV0D7qY3k5ERORYkhDCqr+8P/30E6qrq02uM3jwYGRlZWHhwoW4fPmy/vO2tjYEBARgy5YtuO+++yw6XktLC/r06YN169YhLS0NI0eOxA8//ADpl/LIQghoNBr4+PjghRdewNKlS6/aR3NzM5qbm/Xv6+rqEBUVhdraWgQFBVk0Dmvll1cj7f3vzK63ec44JA7pK8sYiIiIPEldXR2Cg4O79ffb19oNQkNDERoaana9xMRE1NTUoLCwEPHx8QCAvLw8aDQaJCQkWHw8IQSEEPrA5f/+7//Q2NioX37gwAE8/vjj2LNnD4YMGWJwH/7+/vD397f4mPbA9HYiIiLXY3XgY6nhw4cjJSUFc+bMwXvvvYfW1lbMnz8fM2fORGRkJACgsrISSUlJ2LhxI8aOHYtTp07hk08+weTJkxEaGorz589j+fLlCAwMRGpqKgBcFdxcunRJfzylUinX6ViN6e1ERESuR9aUqE2bNiEmJgZJSUlITU3F+PHjsXbtWv3y1tZWlJWVoaGhAQAQEBCAPXv2IDU1FUOHDsWDDz6I3r17Y9++fVdNlHZ15tLbJQARwQEYGx3iyGERERF5Navn+HgCezwjtEROcRWezDpkdPl7s0axYzsREZGF7PH32zWK4BARERE5AAMfmei6tRsjgensREREjsbARybs1k5EROR6GPjIhOnsREREroeBj0yYzk5EROR6GPjIhOnsREREroeBj0xMdWvXvV8yLRY+CmOhEREREdkbAx8ZpcRFIOOOaEhdYhtJAjLuiGYNHyIiIgdj4COjnOIqrN1dga4Z6xoBrN1dgZziKucMjIiIyEsx8JGJro6PqSo9rONDRETkWAx8ZMI6PkRERK6HgY9MWMeHiIjI9TDwkQnr+BAREbkeBj4yMVfHBwBCevohflAfh42JiIjI2zHwkYmpOj466iutmLByJ7O7iIiIHISBj4xS4iKwZtYohAcbf5ylqm3C3KxDDH6IiIgcgIGPzFLiIrDr+YkI6dnD4HJdMjtT24mIiOTHwMcBCs9chvpKi9HlTG0nIiJyDAY+DsDUdiIiItfAwMcBmNpORETkGhj4OIC51HYJQERwAMZGhzhyWERERF6HgY8DmEttFwBS48JRUKHmBGciIiIZSUIIr/tLW1dXh+DgYNTW1iIoKMhhx80prsLS7JJOPbwUEjp1b48IDsCSabFIiYtw2LiIiIjcgT3+fjPwcWDgA2i7thdUqJFbosIHe09ftVx3R2jNrFEMfoiIiDqwx99vPupyMB+FhLHRIdhWrDK4nHV9iIiI5MPAxwkKKtSdHnd1xbo+RERE8mDg4wSs60NEROQcDHycgHV9iIiInIOBjxOYq+sDACE9/RA/qI/DxkREROQNGPg4gbm6PgCgvtKKCSt3sms7ERGRHTHwcZKUuAismTUK4cHGH2epapswN+sQgx8iIiI7YeDjRClxEdj1/ESE9OxhcDlT24mIiOyLgY+TFZ65DPWVFqPLmdpORERkPwx8nIyp7URERI7DwMfJmNpORETkOAx8nMyS1HZloB80QnCeDxERUTfJFvio1Wqkp6cjKCgISqUSs2fPRn19vcltnnjiCQwZMgSBgYEIDQ3F9OnTcezYsavW+/DDD3HzzTcjICAAYWFhmDdvnlynITtLUttrGluRvm4/xq/IY4YXERFRN8gW+KSnp+Po0aPIzc3F1q1bsXv3bmRkZJjcJj4+Hhs2bEBpaSm+/vprCCEwefJktLe369d5/fXX8cILL2Dx4sU4evQotm/fjrvuukuu03AIS1LbAaa3ExERdZckhLD785PS0lLExsbiwIEDGD16NAAgJycHqampOH/+PCIjIy3az5EjRzBixAicPHkSQ4YMweXLl3HttdciOzsbSUlJNo/PHm3t5dCuEfiuvBrzPjqEmsZWg+tIAMKDA/DtoknwUZh6QEZERORZ7PH3W5Y7Pvn5+VAqlfqgBwCSk5OhUCiwf/9+i/Zx5coVbNiwAdHR0YiKigIA5ObmQqPRoLKyEsOHD8eAAQMwY8YMnDt3zuS+mpubUVdX1+nlinwUEhQKyWjQAzC9nYiIqDtkCXxUKhXCwsI6febr64uQkBCoVCqT27777rvo1asXevXqhW3btiE3Nxc9emgL/J06dQoajQZ//etfsXr1anz66adQq9X4zW9+g5YW47Vwli1bhuDgYP1LF0i5Iqa3ExERyceqwGfx4sWQJMnky9BkZGukp6fj8OHD2LVrF4YNG4YZM2agqUn7R16j0aC1tRVvvvkm7rrrLowbNw6bN2/GiRMnsHPnTqP7zMzMRG1trf5l7g6RMzG9nYiISD6+1qy8cOFCPProoybXGTx4MMLDw3Hx4sVOn7e1tUGtViM8PNzk9rq7Mtdffz3GjRuHPn364PPPP0daWhoiIiIAALGxsfr1Q0ND0a9fP5w9e9boPv39/eHv72/m7FyDLr1dVdsEQ5OvdHN8xkaHOHpoREREbs+qwCc0NBShoaFm10tMTERNTQ0KCwsRHx8PAMjLy4NGo0FCQoLFxxNCQAiB5uZmAMBtt90GACgrK8OAAQMAaNPmL126hEGDBllzKi5Ll94+N+sQJKBT8KObyrxkWiwnNhMREdlAljk+w4cPR0pKCubMmYOCggLs3bsX8+fPx8yZM/UZXZWVlYiJiUFBQQEA7fydZcuWobCwEGfPnsW+ffvwwAMPIDAwEKmpqQCAYcOGYfr06fjDH/6Affv2obi4GI888ghiYmIwceJEOU7FKVLiIpBxRzSkLrGNJAEZd0QjJS7COQMjIiJyc7LV8dm0aRNiYmKQlJSE1NRUjB8/HmvXrtUvb21tRVlZGRoaGgAAAQEB2LNnD1JTUzF06FA8+OCD6N27N/bt29dpovTGjRuRkJCAqVOnYsKECfDz80NOTg78/PzkOhWHyymuwtrdFehaqFkjgLW7K1jHh4iIyEay1PFxda5axwfQ1vIZvyIPVbWGs7ZYx4eIiLyVy9bxIdsVVKiNBj0A6/gQERF1BwMfF8M6PkRERPJh4ONiWMeHiIhIPt4d+Fy+7OwRXEVXx8fY7B0JQATr+BAREdnEuwOfUaOAdesAjcbZI9HT1fExNuNcgHV8iIiIbOXdgY9aDcyZAyQkABY2TyUiIiL35d2Bz7JlQFAQcPAgMG4c8NhjwI8/OnVI7RqBpdklRpdLAJZml6C9a5EfIiIiMsu7A5+nngKOHwd0/cc+/BAYNgxYvRpobXXKkJjOTkREJB/vDnwAoH9/YMMGID8fGD0aqKsDFiwAbrkFyMtz+HCYzk5ERCQfBj4648Zp5/m8/z7Qrx9w9CiQlATMmAGY6Pxub0xnJyIikg8Dn44UCuD3v9c+/po/X/t+yxYgJgZ49VWgSf67LExnJyIikg8DH0P69AHeegs4dAi4/XagsRH485+BG28EsrMBGdub6dLZARgMfgSAmWOiZDs+ERGRJ2PgY8qIEcCuXcBHHwGRkcCpU8A99wBTpwInTsh22JS4CKyZNQrhwYYfZ63afgLjV+SxSzsREZGVGPiYI0lAWhpQVgYsXgz4+QHbtgFxcUBmJlBfL8thU+Ii8O2iSViQPMzgclVtE+ZmHWLwQ0REZAUGPpbq1Utb96e4GJgyBWhpAZYv187/2bxZtsdfHx8wPLFadzTW9CEiIrIcAx9rDRsG/PvfwFdfAYMHA5WVwEMPAXfeCRw5YtdDsaYPERGRfTHwsYUkAdOmaVPe//IXIDAQ2L1bW/vn6aft1vyUNX2IiIjsi4FPdwQEAC++CBw7BjzwgLbZ6dtva+8K2aH5KWv6EBER2RcDH3sYOBD417+AHTuA2Fjg0iW7ND81V9MHAJSBftAIwXk+REREFmDgY0+TJgFFRcCqVXZpfmqupg8A1DS2In3dfqa3ExERWYCBj735+QHPPmu35qfmavroML2diIjIPAY+crFj81NdTZ9NsxOgDPQzuA7T24mIiMxj4CM3OzU/9VFIUCgk1DQav2PE9HYiIiLTGPg4gp2anzK9nYiIqHsY+DhSN5ufMr2diIioexj4OIONzU8tTW9va9dg78lL+LKoEvnl1ZzzQ0RE9AtJCJmaTLmwuro6BAcHo7a2FkFBQc4dTH098NprwN//rs346tEDeO454IUXtP3BusgprsLcrEMAfp3QbE5EcACWTItFSlyEHQdORETkWPb4+807Ps5mZfNTS9PbO2KqOxERkRYDH1dhRfNTS9LbO2KqOxERkRYDH1diRfNTS9LbO2KqOxEREQMf12Rh81Nb0taZ6k5ERN6MgY8rM9P81Ja0daa6ExGRN2Pg4w6MND9NePWPiFU0mkxv15Ggze4aGx0i82CJiIhcFwMfd2Gg+aniww/x5VuP47EDX8Kvvc3k5gJAalw4viuvZo0fIiLyWqzj4+w6Prb67jvthOeDBwEA5WHX4cVJc5A/aMRVqyokwFh8wxo/RETkLly6jo9arUZ6ejqCgoKgVCoxe/Zs1NfXm9zmiSeewJAhQxAYGIjQ0FBMnz4dx44d67TOgQMHkJSUBKVSiT59+uCuu+7C999/L9dpuK4uzU+HXDyNzR+/gL2H3sOWKddi0+8T8Pht1wEwHvQArPFDRETeRbbAJz09HUePHkVubi62bt2K3bt3IyMjw+Q28fHx2LBhA0pLS/H1119DCIHJkyejvb0dAFBfX4+UlBQMHDgQ+/fvx7fffovevXvjrrvuQmurZWndHsVA89Nrc7dizJRbkbh5DXYcNt/9nTV+iIjIm8jyqKu0tBSxsbE4cOAARo8eDQDIyclBamoqzp8/j8jISIv2c+TIEYwYMQInT57EkCFDcPDgQYwZMwZnz55FVFQUAOCHH37AzTffjBMnTmDo0KEW7dcjHnUZ8v332sdfe/YAAM4ow/FK0hzsGDJWWyPIjM1zxiFxSF+5R0lERGQTl33UlZ+fD6VSqQ96ACA5ORkKhQL79++3aB9XrlzBhg0bEB0drQ9ybrjhBvTt2xfr169HS0sLGhsbsX79egwfPhzXXXed0X01Nzejrq6u08sjdWh+2hjaH4NqVFj/f3/Bhk9fxnXqSrObs8YPERF5OlkCH5VKhbCwsE6f+fr6IiQkBCqVyuS27777Lnr16oVevXph27ZtyM3NRY8ePQAAvXv3xjfffIOsrCwEBgaiV69eyMnJwbZt2+Dr62t0n8uWLUNwcLD+pQukPJIkAWlpOLJ9P9Yk/BYtCl9MPFWIrz+Yh//Z9SGuaWk0uilr/BARkaezKvBZvHgxJEky+eo6Gdla6enpOHz4MHbt2oVhw4ZhxowZaGrS3olobGzE7Nmzcdttt+G7777D3r17ERcXh6lTp6Kx0fgf9MzMTNTW1upf586d69YY3cHouIHYOP1JpMx+BzsHx8O/vQ1Pffcpdrz/JKaV7Lqq+aky0A9t7ZpOqe4tbRrkl1dbnfrerhE2bUdERCQ3q+b4/PTTT6iurja5zuDBg5GVlYWFCxficofeUm1tbQgICMCWLVtw3333WXS8lpYW9OnTB+vWrUNaWhrWr1+PP/3pT6iqqoJCoei0zvr16zFz5kyL9uuxc3y6yCmuwtysQ4AQmFRegJd2vI9BNdo7bt9FxeHl5CdwLCza6PZd0+AtSX3PKa7C0uwSVNU2WbUdERGROQ6f4xMaGoqYmBiTrx49eiAxMRE1NTUoLCzUb5uXlweNRoOEhASLjyeEgBACzc3NAICGhgYoFApIHSbq6t5rfulfRb9KiYvAmlmjEK4MxI6hCZg8+1387fZZaPT1x7hzxfj3h3/Ay7nvIajJcJmBrjdqzKW+6wKtjkGPJdsRERE5imwFDKdMmYIff/wR7733HlpbW/HYY49h9OjR+OijjwAAlZWVSEpKwsaNGzF27FicOnUKn3zyCSZPnozQ0FCcP38ey5cvx969e1FaWoqwsDAcO3YMI0eOxOOPP46nn34aGo0Gy5cvR3Z2NkpLSxERYdkdBW+546PTrhEoqFDj4s9N6NfTH36V51Az7w+YXKLN/qoODML/TngE/7r5NxCS6VhYAhAeHIBvF02Cj+LXALRdIzB+Rd5VQY+57YiIiCzlslldALBp0ybExMQgKSkJqampGD9+PNauXatf3trairKyMjQ0NAAAAgICsGfPHqSmpmLo0KF48MEH0bt3b+zbt08/UTomJgbZ2dk4cuQIEhMTcfvtt+PChQvIycmxOOjxRj4KCYlD+mL6yGtx2/X90B41EBnTFiFt5ms43ncg+jbWYUXOW/hi40KMvFBmcl8CQFVtEwoq1J0+L6hQGw16TG1HRETkSGxZ4QV3fLr6sqgSf/i4CADg296Ghw/9G89+uwlBLdog9F83JeN/JzyCSz37GN3HGzNHYvrIaw3u05Su2xEREVnKpe/4kOvqmLbe5uOLD8ZMx6SMf2BLXDIAYMYP25G39gk8fuBL+Bppfto19d3SVHimzBMRkTMx8PFCY6NDEBEcgI4zbS717IPnpz6L+2b9DUfChyKopQEv5b2Pf3/4DBLPdO6FZij1fWSUEiE9e5g8rjLQDxohmN5OREROw0ddXvioC+iQ6o5f+3XpKDTtmHEkF8/v3oi+jdoq11tvGI+/TnocF4LCYIipDvBdMb2diIhsYY+/3wx8vDTwAQzX3OlI2VyPZ3dn4XeH/wMfoUGjrz/eSXwA74/9f2j2NX13xxTdnaY1s0Yx+CEiIosx8LERA59fdU11hwRcqm9GWO8AxA/qgwMVarzx+hYs3PoOEs4fBWB981NDmN5ORETWYuBjIwY+lssvr0ba+98BQuCe0t340871CK/XpqTvHByPpUkZOB1ie5YWO8ITEZGlmNVFstN3bJckfBU7AZPm/MPq5qcW7Z+IiMgBGPiQSV3Tzxt6BGLFnY/iLgubn1q7fyIiIjkx8CGTDKW+A0BFyLV47LcvY/b9f8YZZTgi6qvxVvZKfLw5EzEXK8zuV4I2u2tsdIgs4yYiIjKEgQ+Z5KOQsGRaLABcFfxAkow2P13aofmpoanLAsDMMVFyDp2IiOgqnNzMyc0WMZf6rpCA8NqLeCFvPaaW7QUAtPQJwfE/ZCLDdwQu/NxicDvW9CEiIksxq8tGDHxsYy71vfDMZVz8uQlDfyhA7PIXIZWUAADE6NH45OFFWFwZeNU+WdOHiIgsxcDHRgx8HKC1FXjnHWDJEqBOW/3ZWPNT1vQhIiJLMJ2dXJefH/Dss8Dx47h4fxoA481PBYCq2iYUVKidM1YiIvIaDHxIXv37I//FlRY1P2VNHyIikhsDH5JdWO8AHL42Bvf+7u9YfNd8VAcG4YZLZ7H54xfw9hfLEVl3Ub8eERGRnBj4kOx0tYCEwgcfj0zBxIy1+HDU3WiXFLi77FvseH8u/nhgC0RTI1raNMgvr8aXRZXIL69Gu6Ut34mIiCzAyc2c3OwQOcVVmJt1CIB2Tg8ADL94Ci/n/qNT89O/JM3B9g7NT5nuTkREOszqshEDH+cwWAvITPNTprsTEZEOAx8bMfBxnnaNwHfl1Zj30SHUNLbqP7+mpRFP7/sEsw98gR6aNjT7+GL9mHvxduKDaOwRyHR3IiJiOju5Hx+FBIVC6hT0AKabn95dsgtVNY1Mdyciom5j4EMOZypt3VTz06bCQw4cJREReSIGPuRwZtPWjTQ/vTMtBXj6aeDyZccMlIiIPA4DH3I4XXq7udk6zb498PatM5E8Zw3y4u6ApNEAb78NDBsGrFsHaDQOGS8REXkOBj7kcD4KCUumxQKA2eAHACqDwvDDqvfR+vV/0TD0BuDSJWDOHIiEBGD/fnkHS0REHoWBDzlFSlwE1swahfDgzo+9jCVtrdp+AjfsbMHN967AK5PmoK7HNZAOHgTGjQMefxz48UcHjJqIiNwd09mZzu5U7RqBggo1Lv7chLDeAYgf1AdrvinHqu3HTW7X78plLPrmn3igeLv2g6AgYOlSYN48bYNUIiLyOKzjYyMGPq6rXSMwfkVe5yKHJtxSeQx/3bkWwyt/CZRuvBF4801g0iQZR0lERM7AOj7kcQoq1BYHPQBw+NoYTH1oJcpffR3o2xc4ehRISgJmzADOnpVxpERE5I4Y+JBLMVXjxxiNwgfFU2cAJ04A8+cDCgWwZQsQEwO89hrQZP0+iYjIMzHwIZditsaPqe369AHeegs4dAi4/XagsRF48UXt46/sbMD7nuoSEVEXDHzIpVha46ejkJ5+iB/U59cPRowAdu0CPvoIiIwETp0C7rkHmDpVe1eIiIi8FgMfcinW1vgBAPWVVkxYuRM5xVW/fihJQFoaUFYGLFqkzfTatg2IiwMyM4H6evsPnoiIXB4DH3I51tb4AQBVbRPmZh3qHPwAQK9ewPLlQHExMGUK0NKifR8TA3z8MR9/ERF5GaazM53dZXWt8TMySonbVuRBfaXF4PoSgPDgAHy7aBJ8DEVJQgBbtwLPPqt9/AUAEyZo099vvlm28yAiIvtw+XR2tVqN9PR0BAUFQalUYvbs2ai38BGDEAJTpkyBJEn44osvOi07e/Yspk6dimuuuQZhYWF4/vnn0dbWJsMZkDP5KCQkDumL6SOvReKQvig6V2M06AEAAaCqtgkFFWrDK0gSMG2aNuX9L38BAgO1c4FuuYXNT4mIvISsgU96ejqOHj2K3NxcbN26Fbt370ZGRoZF265evRqSdPV/tbe3t2Pq1KloaWnBvn378M9//hMffvghXnrpJXsPn1yMpanuZtcLCNBmex07Bvz2t9pmp2x+SkTkFWQLfEpLS5GTk4N169YhISEB48ePx1tvvYWPP/4YFy5cMLltUVER/v73v+ODDz64atl///tflJSUICsrCyNHjsSUKVPwl7/8Be+88w5aWozfDSD3Z2mqu8Up8QMHauv9bN8OxMbqm5+CzU+JiDyWbIFPfn4+lEolRo8erf8sOTkZCoUC+038UWloaMBDDz2Ed955B+Hh4Qb3e9NNN6F///76z+666y7U1dXh6NGjBvfZ3NyMurq6Ti9yP5akuisD/aARAu0aK6auJSUBRUXA669re36x+SkRkceSLfBRqVQICwvr9Jmvry9CQkKgUqmMbrdgwQLceuutmD59utH9dgx6AOjfG9vvsmXLEBwcrH9FRUVZcyrkIixJda9pbEX6uv0YvyLv6gwvU/z8gAULgOPHgUcf1X62YYP28dfq1UBra3eGTkRELsLqwGfx4sWQJMnk69ixYzYN5quvvkJeXh5Wr15t0/bGZGZmora2Vv86d+6cXfdPjmMs1b0ro+nt5vTvrw148vOB+Higrk4bEN1yC5CX142RExGRK/C1doOFCxfiUd1/ERsxePBghIeH4+LFi50+b2trg1qtNvgICwDy8vJQXl4OpVLZ6fP7778ft99+O7755huEh4ejoKCg0/Iff3kcYWy//v7+8Pf3Nzlmch8pcRH4TWw4viuvxryPDqGm8eq7MQLau0JLs0vwm9hww+ntpowbp53n88EH2oKHuuanDzwA/O1v2vlBRETkdqy+4xMaGoqYmBiTrx49eiAxMRE1NTUoLCzUb5uXlweNRoOEhASD+168eDGOHDmCoqIi/QsAVq1ahQ0bNgAAEhMT8cMPP3QKqnJzcxEUFITY2FhrT4fclI9CgkIhGQx6dMymt5s9iI92sjObnxIReQzZ5vgMHz4cKSkpmDNnDgoKCrB3717Mnz8fM2fORGRkJACgsrISMTEx+js44eHhiIuL6/QCgIEDByI6OhoAMHnyZMTGxuJ3v/sdvv/+e3z99dd48cUXMW/ePN7V8TJ2S283x1Tz061bu7dvIiJyKFnr+GzatAkxMTFISkpCamoqxo8fj7Vr1+qXt7a2oqysDA0NDRbv08fHB1u3boWPjw8SExMxa9YsPPzww3jllVfkOAVyYXZPbzfHUPPTadPY/JSIyI2wZQVbVritdo3A+BV5UNU2wdBFbLaFRXfU1wOvvqpNgW9tBXr0AJ57DnjhBW1/MCIisjuXb1lBJCdz6e0CwMwxMpUu6Nj8NCWFzU+JiNwEAx9ya+bS21dtP2F9TR9rDBsG/Oc/wJdfAoMHA5WVQFoaMHEicOSIPMckIiKbMfAht5cSF4FvF03CguRhBpfbXNPHUpIE3HMPm58SEbkBBj7kMT4+cNbg57qHTkuzS6xrZWEtNj8lInJ5DHzIIxRUqFFVazxtvds1fazB5qdERC6LgQ95BIfV9LEGm58SEbkcBj7kERxe08dSbH5KRORSGPiQRxgbHYKI4ACjXdsBQBnoB40Q8s7zMYbNT4mIXAIDH/II5mr6AEBNYyvS1+2XN73dHF3z07Vrgb59f21+OmMGcNbw5GwiIrIfBj7kMczV9NGRPb3dHF3z0+PHgXnz2PyUiMiBGPiQR9HV9Nk0OwHKQD+D6zgsvd2ckBBtujubnxIROQwDH/I4PgoJCoWEmkbjE4cdmt5uDpufEhE5DAMf8kgumd5uiiRpW12UlQGLFmmzwf7zHyAuDsjM1DZFJSKibmPgQx7JZdPbzWHzUyIiWTHwIY9kLr1dAhARHICx0SGOHJbl2PyUiEgWDHzII5lKb9e9XzItFj4KU5V/nIzNT4mI7I6BD3mslLgIZNwRDalLbCNJQMYd0UiJi3DOwKzF5qdERHbDwIc8Vk5xFdburkDXjHWNANburnBeHR9bsfkpEVG3MfAhj9SuEViaXQJTU4GdXsfHVmx+SkRkMwY+5JEKKtSoqjWequ5SdXxsoWt+WlZ2dfPTN95g81MiIiMY+JBHcrs6PrYKD7+6+emzz2onQO/c6ezRERG5HAY+5JHcto6PrQw1P500ic1PiYi6YOBDHslcHR8AUAb6oa1dg70nL+HLokrkl1e755wfHTY/JSIySxLC+0rB1tXVITg4GLW1tQgKCnL2cEgmOcVVmJt1CABMTnLuKCI4AEumxbpPqrsp33+vrfezZ4/2/eDB2vk/d9/t3HEREdnIHn+/eceHPFZKXATWzBqF8GDLH2epapswN+uQ+6W6G8Lmp0REV2HgQx4tJS4C3y6ahE2zE6AM9DO7vu7OkNumunfF5qdERJ0w8CGP56OQoFBIqGm0LMXb7VPdDWHzUyIiAAx8yEvYkrbu9qnuhnRsfhodzeanROR1GPiQV7Albd1jUt270jU/LSm5uvnpM8+w+SkReTQGPuQVLElv70gZ6AeNEJ4xz8cYQ81P33qLzU+JyKMx8CGv4KOQsGRaLABYFPzUNLYifd1+jF+R5xkZXqaw+SkReREGPuQ1vD693Rw2PyUiL8AChixg6HXaNQIFFWpc/LkJ/Xr6QyMEnt582GjWlwQgPDgA3y6aBB+FpQ/L3JxKpU13//BD7fugIOCVV4CnntKmxBMROQELGBLZwEchIXFIX0wfeS1uu74ffH0UJlPdPTK93Rw2PyUiD8XAh7ye13RytwWbnxKRh2HgQ17P6zq5W4vNT4nIg8ga+KjVaqSnpyMoKAhKpRKzZ89GvYUl8oUQmDJlCiRJwhdffKH//Pvvv0daWhqioqIQGBiI4cOH44033pDpDMgbmEt1l6BtXjo2OsSRw3I9ISHA228Dhw4Bt98ONDZq0+FvvBHYutXZoyMisoisgU96ejqOHj2K3NxcbN26Fbt370ZGRoZF265evRqSdPWfosLCQoSFhSErKwtHjx7FCy+8gMzMTLz99tv2Hj55CXOp7gLAzDFRDh2TS9M1P920ic1PicjtyJbVVVpaitjYWBw4cACjR48GAOTk5CA1NRXnz59HZGSk0W2Liopw99134+DBg4iIiMDnn3+Oe++91+j68+bNQ2lpKfLy8iwaG7O6yJCc4ioszS5BVa3hRzcRwQFYMi0WKXERDh6ZC6uvB159VZsC39oK9OgBPPcc8MIL2v5gRER25NJZXfn5+VAqlfqgBwCSk5OhUCiw30RRtIaGBjz00EN45513EB4ebtGxamtrERJi/DFEc3Mz6urqOr2IutJ1cl+QPMzgcq+q6WMpNj8lIjcjW+CjUqkQFhbW6TNfX1+EhIRApVIZ3W7BggW49dZbMX36dIuOs2/fPnzyyScmH6EtW7YMwcHB+ldUFB9bkHEfHzCcraT7E740u8SzW1nYgs1PichNWB34LF68GJIkmXwdO3bMpsF89dVXyMvLw+rVqy1av7i4GNOnT8eSJUswefJko+tlZmaitrZW/zp37pxN4yPPV1ChNvqoC/DSmj6WYvNTInIDvtZusHDhQjz66KMm1xk8eDDCw8Nx8eLFTp+3tbVBrVYbfYSVl5eH8vJyKJXKTp/ff//9uP322/HNN9/oPyspKUFSUhIyMjLw4osvmhyPv78//P39Ta5DBLCmj13omp8+/DCwcCHw6afa5qebNwPLlmlbYChYSYOInMPqwCc0NBShoaFm10tMTERNTQ0KCwsRHx8PQBvYaDQaJCQkGNxm8eLF+P3vf9/ps5tuugmrVq3CtGnT9J8dPXoUkyZNwiOPPILXXnvN2lMgMoo1fexI1/x0xw7tHZ+SEm09oH/8Q5sWb+TfASIiOcn2n13Dhw9HSkoK5syZg4KCAuzduxfz58/HzJkz9RldlZWViImJQUFBAQAgPDwccXFxnV4AMHDgQERHRwPQPt6aOHEiJk+ejOeeew4qlQoqlQo//fSTXKdCXsRcTR8AUAb6oa1dg70nL+HLokrkl1ejpU2D/PJqo++9ek4Qm58SkQuRtUmpWq3G/PnzkZ2dDYVCgfvvvx9vvvkmev2S5nr69GlER0dj586duPPOOw0PUJI6pbO//PLLWLp06VXrDRo0CKdPn7ZoXExnJ1NyiqswN+sQgF8nNJujkICOsU3X90yF/wWbnxJRN9jj7ze7szPwIQPM1fSxlu4O0ppZoxj8ANrmp08/DRQWat/feKN2HtDEic4dFxG5NJeu40PkznQ1fTbNToAysPt3IpgK30ViIpufEpFTMPAhMsJHIUGhkFDT2GqX/TEVvgs2PyUiJ2DgQ2SCHGnrTIXvgs1PiciBGPgQmSBH2jpT4Y1g81MicgAGPkQmWJLebikJ2uyusdHG+8p5PUkCHnoIKCsDFi3SZnr95z9AXJw2G6y+3tkjJCI3x8CHyAQfhYQl02IBoNvBjwAwc0wU2jWCNX7MYfNTIpIJ09mZzk4WMJfebq6Oj6l1WePHDCGA7Gzg2WeBigrtZxMmAG++Cdx8s1OHRkSOxTo+NmLgQ7Zo1wgUVKhx8ecm9OvpD0jApfpmhPUOQPygPig8cxkXf27Sv1/zTTlWbT9udr+s8WOhpiZg5Uptv6/GRm0W2Lx5wNKlQJ8+zh4dETkAAx8bMfAhubVrBMavyLO4AKIEIDw4AN8umgQfhT1mFHmwM2eAP/5R2/wUAPr1Y/NTIi/BAoZELqqgQm1V1WfW+LHCoEHaej/btwOxscClS9p6QAkJ2qKIREQmMPAhkoGttXpY48cKbH5KRDZg4EMkA1tr9bDGj5X8/IAFC7Tp748+qv1swwZg2DDgjTeAVvtU3SYiz8HAh0gGttT/Cenph/hBnKRrk/BwbcCzbx8QHw/U1WmzwG65Bdi509mjIyIXwsCHSAa21P9RX2nFhJU7kVNcJd/APB2bnxKRGQx8iGSSEheBNbNGITy48+MrU0lbqtomzM06xOCnO9j8lIhMYDo709lJZh3r/4T1DsDIKCVuW5EH9ZUWg+sztd3Ovv8eePppYM8e7fvBg7Xzf+6+27njIiKrMZ2dyA34KCQkDumL6SOvReKQvig6V2M06AGY2m53bH5KRB0w8CFyMEtT1pnabke65qfHjrH5KZGXY+BD5GCWpqwztV0GvXuz+SmRl2PgQ+RglqS6KwP90Nauwd6Tl/Rd3FvaNJ26und9zy7vVhg2THvH58svgehooLISSEsDJk4Ejhxx9uiISEac3MzJzeQEOcVVmJt1CIB2To8lzHWAZ5d3G7H5KZHb4ORmIjdlLNXdlK43dLq+Zyq8jQICgD//GSgtBX77W0CjAd56S3tXaN067Xsi8hgMfIicJCUuAt8umoRNsxOgDPTr9v50cdDS7BI+9rKFsean48ax+SmRB2HgQ+REPgoJCoWEmkb79JRiKrwddG1+euAAm58SeRAGPkROJkfaOlPhu4nNT4k8FgMfIieTI22dqfB2wuanRB6HgQ+Rk9nSyd0YCdrsrrHRIXbYG+mx+SmRx2A6O9PZyQXYkt5uzILk6zH3zqEoPHMZF39uQr+e/oAEXKpvRljvAMQP6tPtZYbej40O8Y7eYmo18NJLwJo12oyvwEDghReAhQu1GWJEJBt7/P1m4MPAh1xETnEVlmaXoKrW8Pwcc3V8TK0rxzKvryPE5qdEDsfAx0YMfMhVdezkbskdlzXflGPV9uPOHjYA6B/VrZk1ynuCHyGAzZuB558HLlzQfpaaCqxeDVx/vVOHRuSJGPjYiIEPeYJ2jcD4FXlG7xA5gwQgPDgA3y6a5B2PvXR+/hl47TVtCnxrK9Cjh/bR15/+BPTq5ezREXkMVm4m8mIFFWqXCnoAL64jZKj56bJlbH5K5IIY+BC5KVeu1ePKY5MVm58SuTwGPkRuypVr9bjy2GQnScA99wAlJcArr2izvnbt0tb+eeYZ4PJlZ4+QyKsx8CFyU/as/2NPIT39MDJKifzyanxZVIn88mrv7B3G5qdELknWwEetViM9PR1BQUFQKpWYPXs26uvrLdpWCIEpU6ZAkiR88cUXBteprq7GgAEDIEkSampq7DdwIjfgo5CwZFosALhU8KO+0oobl+Qg7f3v8IePi5D2/ncYvyLPe7vGd2x+Onw4m58SOZmsgU96ejqOHj2K3NxcbN26Fbt370ZGRoZF265evRqSZPqf89mzZ+Pmm2+2x1CJ3FJKXATWzBqF8ODOj5ZMJVTZa5mpdbve4FHVNmFu1iHvDX4AbfPT779n81MiJ5Mtnb20tBSxsbE4cOAARo8eDQDIyclBamoqzp8/j8jISKPbFhUV4e6778bBgwcRERGBzz//HPfee2+nddasWYNPPvkEL730EpKSknD58mUolUqLxsZ0dvI0Hev/2Ks6s7k6QiOjlLhtRR7UV1osGqPXprobolIBmZnAhx9q3wcFaecDPfWUtkEqERlkj7/fvnYek15+fj6USqU+6AGA5ORkKBQK7N+/H/fdd5/B7RoaGvDQQw/hnXfeQXh4uMF1SkpK8Morr2D//v04deqU2bE0NzejublZ/76urs7KsyFybT4KCYlD+nb6rOt7ey/LL6+2OOgBOqe6mzqGV9A1P83I0FZ/LizUNj99/33tPKCJE509QiKPJdujLpVKhbCwsE6f+fr6IiQkBCqVyuh2CxYswK233orp06cbXN7c3Iy0tDSsXLkSAwcOtGgsy5YtQ3BwsP4VFRVl+YkQkUG2pqx7baq7Icaanz74IHDunLNHR+SRrA58Fi9eDEmSTL6OHTtm02C++uor5OXlYfXq1UbXyczMxPDhwzFr1iyL95uZmYna2lr96xz/QSHqNltT1r061d0QHx/tZOfjx4F58wCFAvjXv7TFD197DWhioEhkT1YHPgsXLkRpaanJ1+DBgxEeHo6LFy922ratrQ1qtdroI6y8vDyUl5dDqVTC19cXvr7aJ3H3338/7rzzTv06W7Zs0S9PSkoCAPTr1w9LliwxuF9/f38EBQV1ehFR91ibTi9B28h0bHSInMNyXyEhwNtvA4cOAbffDjQ0AC++CNx4I7B1q7NHR+QxZJ/cfPDgQcTHxwMA/vvf/yIlJcXo5GaVSoVLly51+uymm27CG2+8gWnTpiE6Ohrl5eVobGzULz9w4AAef/xx7Nu3D0OGDLnq8ZohnNxMZB85xVWYm3UIgHYOjzkLkq/H/EnXc3KzOWx+SmSQyzcpnTJlCn788Ue89957aG1txWOPPYbRo0fjo48+AgBUVlYiKSkJGzduxNixYw0PUJIMZnXpfPPNN5g4cSKzuoicJKe4CkuzSzr1DVNIV6e060QEB2DJtFjv6eDeHWx+StSJyzcp3bRpE2JiYpCUlITU1FSMHz8ea9eu1S9vbW1FWVkZGhoa5BwGEckoJS4C3y6ahM1zxuGNmSOxec44HPvLFCxIHmZwfdb0sYKu+ekPP7D5KZGdyHrHx1Xxjg+RvNo1AuNX5BntHs+aPjYQAsjO1qa9V1RoP5swAXjzTYCFXMlLuPwdHyLyTgUVaqNBD9C5pg9ZiM1PieyCgQ8R2Z2ltXpY08cGbH5K1C0MfIjI7iyt1cOaPt1gqvlpQYGzR0fkshj4EJHdWVLjRxnoB40QaDeW/kWWMdT8NCGBzU+JjGDgQ0R256OQsGRaLAAYDX5qGluRvm4/xq/IY4ZXd/n5AQsWAGVlwKOPaj/bsEH7+OuNN7Sp8EQEgIEPEckkJS4Ca2aNQniw6cdZTG+3I13z0337gPh4oK5OmwV2yy3Azp3OHh2RS2DgQ0Sy0dX42TQ7AcpAP4Pr6B50Lc0u4WMve9E1P/3HP9j8lKgLBj5EJCsfhQSFQkJNo/HHLUxvl4GPD5CRweanRF0w8CEi2TG93YnY/JSoEwY+RCQ7pre7gBEjtAUPN20CIiOBU6eAadOAqVOBEyecPToih2HgQ0SysyS9PaSnH+IH9XHYmLySJAEPPQQcOwYsWqTNBvvPf4C4OG3j0/p6Z4/QpbRrBPLLq/FlUSXyy6vR0qbRv9974hL2nrxkcFl+eTXnq7kw9upiry4ih8gprsLcrEMAfp3Q3BU7tztYWZk26ysnR/v+2muBv/1NOwla8u4eajnFVViaXdKp9YpCAozFM12X8VqWB3t1EZHbsCS9nantDnbDDdo7Pl9+CURHA5WVQFoaMHEicOSIs0fnNLogvWu/OVM3cbou47Xsuhj4EJHDpMRFYNfzExHSs4fB5UxtdwI2P+2kXSOwNLvE6F1JS/Fadl0MfIjIoQrPXIb6SovR5UxtdxI2PwUAFFSor7rTYytey66JgQ8RORRT212clzc/leO647XsWhj4EJFDMbXdTXhp81M5rjtey66FgQ8ROZS51HYJ2oyYsdEhjhwWGeKFzU8tKb1gKV7LromBDxE5lLnO7QJAalw4CirUnBTqKjyo+amp2jy6+jszxwzs9uRmQHstzxwTZdExLRkrfx/sg3V8WMeHyCksqZPCWiguqL0dWL9eW/Cwulr72YwZ2vo/UVHOHZsZllxz1tTqsXWZJde5obHy98E+f78Z+DDwIXKado1AQYUauSUqfLD39FXLdXeE1swa5dX/2LsktRp46SVgzRptxtc112iDoYULtRliLkZXm8fWP3gLkq/H3DuHovDMZVz8uQn9evoDEnCpvhlhvQMQP6gP1nxTjlXbj1u9767XubGx8veBgY/NGPgQuY52jcD4FXlGU4glAOHBAfh20ST4KLy7mrBL+v574OmngT17tO8HD9bO/7n7bueOqwNz15g5llyD9jrGrucnYsLKnfx9MIKVm4nI7Zmrm8JaKC7ODZqfdrc2jyXXoL2O8f/ln+bvg8wY+BCRU7Gujwdw8ean9rp2TO3HXsc4o27o9ljINAY+RORUrOvjQXr3BpYvB374AUhJAVpagGXLgJgY4OOPASfNrLDXtWNqP/Y6xqCQa7o9FjLN19kDICLvpquboqptMjrxVBnoh7Z2DfaevKSfTDo2OsQj5jjoJnhf/LnJc85L1/w0O1ub9l5RoW1++t57wJtvAjffbNNuO/6sDE0uNjbxeGSUEiE9e5hslWKKbl6NqXo8llzH5oT09MODYwbinW/KTY5VGegHjRBo1wj3v1acgJObObmZyOl0WSwALP6j4QmpvV6RstzUBKxcqb3z09gIKBTAvHnA0qVAnz4W78bQz6ojW1PNzbEmk8qW67gra8bqcdeKBTi5mYg8QkpcBNbMGoXwYMtv36tqmzA36xByiqtkHJl8dH8ku/4hd/fzuoodmp8a+1l1ZCpYMLWs6w2Tru/DgwMsTh83dh2bO4alY+3K464VB+EdH97xIXIZ7RqB78qrMe+jQ6hpNN8OwV1Te706hX/HDm36e2mp9v2YMcDbbwNjxxpcvbtp4qaE9PTD3kVJKDpXo3/U2PGRma2PHrs+vuy6z5FRSty2Is/mR28defS1YoA9/n5zjg8RuQwfhQSFQrIo6AE6p/YmDukr7+DsyJoUfnc6L4vomp++/Tbw8su/Nj997DHt47D+/Tut3t00cVPUV1pRdK7mqp9xd3/mPgrJ5D7zy6vtEvQAHn6tyISPuojIpdiSputuqb1en8JvRfNTuX8GzvgZy3FMj71WZMDAh4hcii1puu6W2ssU/l9Y0PxU7p+BM37GchzT468VO+KjLiJyKdamBRtKdbfHPA1zLE2tNpRqrRECykA/k4/0Qnr6YWSUEvnl1QaP4Wpp75am5Ruc/xI2DBf/8Rlu/M+nGLTqNfgdPQpMmoT2Bx6A7/wXzP6sbGFJirpc7JH63lHXa8Ve14e5uUqmSgh0PL6pa8OaY4T1DkBM3+6HLZzczMnNRC6nu2nBcnd5tza12tZ0alPbuVIqs6Vp+ZZ0Rw9u/BnPfZuFWYe3wUdo0ODnj3fGzcC6sfeh2beHXcbrCs0+Lb3GpQ7LJRPr2vua724ne93xARi9NgwtM3eMUP92HHxlOpuUWouBD5HrMxdcWMOef+i62+XbXlzhjzdg/OdhacdxY4ZfPIWXc/+BhPNHAQCnlRF4JWkO8oYazv6ypo6PqwSNlgQXpoIEU7pzfdjjGjcVpJlaZo5obsDZ1TNcM/BRq9V4+umnkZ2dDYVCgfvvvx9vvPEGevXqZXZbIQRSU1ORk5ODzz//HPfee2+n5R9++CFef/11HD9+HEFBQXjggQfwzjvvWDw2Bj5E7qHr4ySNEHh682GbHnvYI+1XztRqWzg7ldnStHxzHceNEgL3lO7Cn3Z+gPB6bVPOvMGjsXrqU3h+/t1QKCSLKjc74tGnrcw96uk41pY2DcYt22FxRpgt14erXeNdaZobcK6bgY9sc3zS09NRVVWF3NxctLa24rHHHkNGRgY++ugjs9uuXr0akmT4S3r99dfx97//HStXrkRCQgKuXLmC06dP23n0ROQKuqYF55dX2zzXwx5pv3KmVtvC2anMlqblm+s4bpQk4avYO7FjyFjMz/8XZh/4ApNOHcRt72bgp9b5GLDyL0CvfvrVTf0MXDXV21zqe0eFZy5blQZvy/Xhate4HGQJfEpLS5GTk4MDBw5g9OjRAIC33noLqamp+Nvf/obIyEij2xYVFeHvf/87Dh48iIiIzrfnLl++jBdffBHZ2dlISkrSf36zjX1fiMi92CNltzv7cNWUYWeNy9LjWtpx3Jgr/tdgxZ2PYstNyXhpx/u4s6IQA95bDWRvAf72N+DBB7Ud4j2crd+zNdu56jVuT7Kks+fn50OpVOqDHgBITk6GQqHA/v37jW7X0NCAhx56CO+88w7Cw8OvWp6bmwuNRoPKykoMHz4cAwYMwIwZM3Du3DmT42lubkZdXV2nFxG5H3uk7HZnH66aMuyscVl6XEs7jptzqu8APPrAy/j9//szmqIGAZWV2uanEycCR47Y5RiuzNbv2ZrtXPUatydZAh+VSoWwsLBOn/n6+iIkJAQqlcrodgsWLMCtt96K6dOnG1x+6tQpaDQa/PWvf8Xq1avx6aefQq1W4ze/+Q1aWozf/lu2bBmCg4P1r6ioKNtOjIicSpcGbMt/20vQThTtTvpyd44vB3ucU3eY+3noxve7xOvs9nOTJAlHR0+A37FS4JVXgMBAYNcube2fZ54BLl+2w1Fck7XXny3Xh6td43KwKvBZvHgxJEky+Tp27JhNA/nqq6+Ql5eH1atXG11Ho9GgtbUVb775Ju666y6MGzcOmzdvxokTJ7Dzl2JXhmRmZqK2tlb/MneHiIhck49C0me4WPsPswAwc0wU2jUC+eXV+LKoEvnl1Whp0+jf7z1xCXtPXjK67LtT1ZgSFw5hw/GNkYz8f0vPKTUuHAUVaovPI7+8Gu22tirvQDcp19TPQ/czB4CZYwZ2OxNOd4wl02Lhc01gt5ufuhtrr3/d9fFdebXJ67rjsnaNsMt35cqsyur66aefUF1dbXKdwYMHIysrCwsXLsTlDpF3W1sbAgICsGXLFtx3331Xbffss8/izTffhELxayzW3t4OhUKB22+/Hd988w02bNiAxx9/HOfOncOAAQP06/Xv3x+vvvoq5syZY9F5MKuLyL1ZW0fH1DJr1pXjGNbUNLHXeTiixoup49u6zOS4rWx+6s66W2PHEde1XMewR1aXLOnspaWliI2NxcGDBxEfHw8A+O9//4uUlBScP3/e4ORmlUqFS5cudfrspptuwhtvvIFp06YhOjoax48fxw033IDt27frJzer1WqEhoZi27ZtmDx5skXjY+BD5P7MVU5e8005Vm0/Lvs4Zt92HSbF9LeqcrOlFW4Nbbej9Ed8sPd0t8bs7BovOguSr8fcO4faLw29tVUb7CxZAvz8s/azxx/XNj/tMv3C3Rm7VnJLVN2+Pgzp+l3JcX0aO0bH32WXDXwAYMqUKfjxxx/x3nvv6dPZR48erU9nr6ysRFJSEjZu3IixRiJySZKuquNz77334uTJk1i7di2CgoKQmZmJU6dOoaioCH5+fhaNjYEPkWdzVC0SR9fRsed5ObvGi6w/O5UKWLwY+Oc/te+Dg4GlS4GnntI2SPVQcl33ln5X3Tm+sWN03ac9Ah/ZmpRu2rQJMTExSEpKQmpqKsaPH4+1a9fql7e2tqKsrAwNDdalOW7cuBEJCQmYOnUqJkyYAD8/P+Tk5Fgc9BCR53NULZKOdVIcwZ7nZcvYnX18i4WHAx9++Gvz09raq5qfeiK5rntLv6vuHN/YMeQ4J9kKGIaEhJgsVnjdddfB3M0mQ8uDgoKwfv16rF+/vttjJCLP5OhaJI46nhzHcXaNF1l/domJwP79wPr1wJ/+BPzS/BQzZmjr/3hYhq/c16G5/ctRZ0uOc5Ltjg8RkbM4uhaJo44nx3GcXeNF9p+djw+QkQEcPw7MmwcoFMC//gXExACvvQY0eU7BPrl/lub2L0edLTnOiYEPEXkcR9UicXQdHXufV0hPP4yMUlqcBj8ySomQnvbrkO7QGkQhIdqJz4cOAbffDjQ0AC++CMTFAVu3OmYMMpPrurf0u5KjzpYc58Tu7JzcTOSRdNlHgO2doE1xVnd0e5+XvdKXO3bcNtd92+md5YUANm8Gnn8euHBB+1lqKrB6NXD99Y4fjx3Z+/qw9rsydnxT14e5Y+QUV+HJX/bp0pObiYicKSUuAmtmjUJ4cOdb5V2TUkwlFJlaFh4c4JQ/3Jael6XLTNUy7LrM1LrhwQF4b9YovGfB2Jz1s9OTJOChh4Bjx4BFi7SZXv/5j/buz5/+BNTXO2dcdmDs+ujImmvF2u/K2PFNXR+Ovh54x4d3fIg8mqnaOOZq7Fhaf8fVzqvrWEdGKXHbijyrOntbI6SnH77LTEYPX4XZsbnCz+4qZWXarK+cHO37a691++an5upc2a12kgXHN1evytQx5Ehn98rAp7a2FkqlEufOnWPgQ0Qer+CUGo//84Csx/jgkTEYO9g5PcPsQghg2zZt/Z8zZ7Sf3XYb8L//q70TRE7R9drVNDegcs2jqKmpQXBwsE379MrA5/z582xUSkRE5Ka6tq6yhlcGPhqNBhcuXEDv3r0hucBtzLq6OkRFRXnsHSien3vj+bk3np/78uRzA2w7PyEEfv75Z0RGRnbq7WkN2QoYujKFQmFzpCinoKAgj7y4dXh+7o3n5954fu7Lk88NsP78bH3EpcOsLiIiIvIaDHyIiIjIazDwcQH+/v5YsmQJ/P39nT0UWfD83BvPz73x/NyXJ58b4Lzz88rJzUREROSdeMeHiIiIvAYDHyIiIvIaDHyIiIjIazDwISIiIq/BwKebli1bhjFjxqB3794ICwvDvffei7Kysk7r3HnnnZAkqdPrySeftPgYTz75JCRJwurVqzt9rlarkZ6ejqCgICiVSsyePRv1du4qLNf5vfzyy4iJiUHPnj3Rp08fJCcnY//+/Z3Wue66667a7/Llyz3i3Nz1u2ttbcWiRYtw0003oWfPnoiMjMTDDz+MCxcudFpP7u/O2efnrt8fAHz22WeYPHky+vbtC0mSUFRUdNU63f03yxLOPL+mpibMmzcPffv2Ra9evXD//ffjxx9/tOfpyXZ+Qgi89NJLiIiIQGBgIJKTk3HixIlO67jr75+l59ft3z9B3XLXXXeJDRs2iOLiYlFUVCRSU1PFwIEDRX19vX6dCRMmiDlz5oiqqir9q7a21qL9f/bZZ2LEiBEiMjJSrFq1qtOylJQUMWLECPHdd9+JPXv2iKFDh4q0tDR7np5s57dp0yaRm5srysvLRXFxsZg9e7YICgoSFy9e1K8zaNAg8corr3Tab8fjuvO5uet3V1NTI5KTk8Unn3wijh07JvLz88XYsWNFfHx8p/Xk/u6cfX7u+v0JIcTGjRvF0qVLxfvvvy8AiMOHD1+1Tnf+zXKH83vyySdFVFSU2LFjhzh48KAYN26cuPXWW93i/JYvXy6Cg4PFF198Ib7//ntxzz33iOjoaNHY2Khfx11//yw9v+7+/jHwsbOLFy8KAGLXrl36zyZMmCD+8Ic/WL2v8+fPi2uvvVYUFxeLQYMGdQp8SkpKBABx4MAB/Wfbtm0TkiSJysrK7pyCSfY8v45qa2sFALF9+3b9Z13PWW6OOjdP++4KCgoEAHHmzBn9Z47+7oRw3Pl5yvdXUVFhMvDp7s/NWo46v5qaGuHn5ye2bNmi/6y0tFQAEPn5+TYdyxL2OD+NRiPCw8PFypUr9Z/V1NQIf39/sXnzZv1n7vr7Z8n52eP3j4+67Ky2thYAEBIS0unzTZs2oV+/foiLi0NmZiYaGhpM7kej0eB3v/sdnn/+edx4441XLc/Pz4dSqcTo0aP1nyUnJ0OhUFz1WMWe7HV+HbW0tGDt2rUIDg7GiBEjOi1bvnw5+vbti1tuuQUrV65EW1tb90/CCEedmyd9d7r9SpIEpVLZ6XNHfne6cQDyn5+nfX/GyLVfYxx1foWFhWhtbUVycrL+s5iYGAwcOBD5+fnd2rcp9ji/iooKqFSqTmMPDg5GQkLCVWN3x98/S87PHr9/XtmkVC4ajQbPPvssbrvtNsTFxek/f+ihhzBo0CBERkbiyJEjWLRoEcrKyvDZZ58Z3deKFSvg6+uLZ555xuBylUqFsLCwTp/5+voiJCQEKpXKPifUhT3PDwC2bt2KmTNnoqGhAREREcjNzUW/fv30y5955hmMGjUKISEh2LdvHzIzM1FVVYXXX3/drc/NE747naamJixatAhpaWmdmgw68rsDHHt+nvT9GSPXfo1x5PmpVCr06NHjqkC9f//+Lv/96cbXv39/k2N3198/S87PLr9/Ft+DIrOefPJJMWjQIHHu3DmT6+3YsUMAECdPnjS4/ODBg6J///6dbtt1vXX52muviWHDhl21bWhoqHj33XdtOwEz7HV+OvX19eLEiRMiPz9fPP744+K6664TP/74o9H1169fL3x9fUVTU5NN4zfFkefmCd+dEEK0tLSIadOmiVtuucXsc3s5vzshHHt+nvL9mXrU1Z392sKR57dp0ybRo0ePq9YfM2aM+J//+R+rxm0pe53f3r17BQBx4cKFTp8/8MADYsaMGUb36y6/f5acnz1+//ioy07mz5+PrVu3YufOnRgwYIDJdRMSEgAAJ0+eNLh8z549uHjxIgYOHAhfX1/4+vrizJkzWLhwIa677joAQHh4OC5evNhpu7a2NqjVaoSHh3f/hLqw5/np9OzZE0OHDsW4ceOwfv16+Pr6Yv369Sb329bWhtOnT1s9flMcfW6e8N21trZixowZOHPmDHJzczvd7TG2Xzm+O8Dx5+cJ35+15Nov4PjzCw8PR0tLC2pqajp9/uOPP7r896cbX9cMNHNjd5ffP0vOzx6/fwx8ukkIgfnz5+Pzzz9HXl4eoqOjzW6jS6+MiIgwuPx3v/sdjhw5gqKiIv0rMjISzz//PL7++msAQGJiImpqalBYWKjfLi8vDxqNRn9x2YMc52eMRqNBc3Ozyf0qFIqrbnPaylnn5u7fnS4oOHHiBLZv346+fftatF97fneA887P3b8/W8ixX2edX3x8PPz8/LBjxw79Z2VlZTh79iwSExNt3m9XcpxfdHQ0wsPDO429rq4O+/fvNzl2d/n9s+T87PL7Z9F9ITJq7ty5Ijg4WHzzzTedUvYaGhqEEEKcPHlSvPLKK+LgwYOioqJCfPnll2Lw4MHijjvu6LSfG264QXz22WdGj2Noln5KSoq45ZZbxP79+8W3334rrr/+erun1MpxfvX19SIzM1Pk5+eL06dPi4MHD4rHHntM+Pv7i+LiYiGEEPv27ROrVq0SRUVFory8XGRlZYnQ0FDx8MMPu/25CeG+311LS4u45557xIABA0RRUVGn/TY3NwshHPPdOfP8hHDf708IIaqrq8Xhw4fFv//9bwFAfPzxx+Lw4cOiqqrKqv266/kJoX00M3DgQJGXlycOHjwoEhMTRWJioluc3/Lly4VSqRRffvmlOHLkiJg+fXqndG93/v2z5PyE6P7vHwOfbgJg8LVhwwYhhBBnz54Vd9xxhwgJCRH+/v5i6NCh4vnnn79qTkTHbQwxFPhUV1eLtLQ00atXLxEUFCQee+wx8fPPP7v8+TU2Nor77rtPREZGih49eoiIiAhxzz33iIKCAv36hYWFIiEhQQQHB4uAgAAxfPhw8de//tWuz6iddW5CuO93p5s3Yei1c+dOIYRjvjtnnp8Q7vv9CSHEhg0bDO53yZIlVu3XXc9PCO3v6VNPPSX69OkjrrnmGnHfffd1Coxc+fw0Go3485//LPr37y/8/f1FUlKSKCsr0y93598/S85PiO7//km/HJiIiIjI43GODxEREXkNBj5ERETkNRj4EBERkddg4ENEREReg4EPEREReQ0GPkREROQ1GPgQERGR12DgQ0RERF6DgQ8RERF5DQY+RERE5DUY+BAREZHXYOBDREREXuP/B21ieCBA7egCAAAAAElFTkSuQmCC",
      "text/plain": [
       "<Figure size 640x480 with 1 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "#Testing the derived skeleton positions\n",
    "Name01 = 'Skeletons/'+List_Names[1]+'_skeleton.fits'\n",
    "Test01 = fits.open(Name01)\n",
    "\n",
    "# Running the function\n",
    "TestPos01 = OF.PositionList(Test01)\n",
    "# Testing the linear fit to the data\n",
    "TestLine01 = OF.LinearFit(TestPos01)\n",
    "#print(TestLine01[1])\n",
    "# Creating predicted Y positions\n",
    "TestPred01 = np.copy(TestPos01)\n",
    "TestPred01[1,:] = TestLine01[0]*TestPred01[0,:] + TestLine01[1]\n",
    "#print(TestPred01[1,:].argmax())\n",
    "\n",
    "# Testing the angle \n",
    "TestAngle01 = OF.CalculateAngle(TestPos01,TestLine01)\n",
    "print(TestAngle01)\n",
    "\n",
    "# Plotting astronomical picture and positions\n",
    "Test01_Plot01 = FITSFigure(Test01)\n",
    "Test01_Plot01.show_colorscale(vmin=0.0,vmax=0.5)\n",
    "Test01_Plot01.show_markers(TestPos01[0,:],TestPos01[1,:],edgecolor='black',facecolor='dodgerblue',marker='*',s=50,alpha=0.8, zorder=2,linewidth=0.4)\n",
    "\n",
    "# Plotting data and fit\n",
    "Test01_Plot02 = plt.figure()\n",
    "Test01_Plot02 = plt.scatter(TestPos01[0,:],TestPos01[1,:])\n",
    "plt.axis([max(TestPos01[0,:]), min(TestPos01[0,:]), min(TestPos01[1,:]), max(TestPos01[1,:])])\n",
    "plt.plot(TestPred01[0,:],TestPred01[1,:], color = 'r')"
   ]
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "PoLiteWIP",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.10.8"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}
